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Ricci scalar is the key ingredient of non-Newtonian theory of grav- 
ity, where space-time geometry has a crucial role. Normally, it is sup- 
posed to be a geometrical field, but interestingly it also behaves like 
a physical field. Thus it plays dual role in the arena of gravitation. 
This article is an overview of the work related to dual roles of the 
Ricci scalar. A scalar is a mathematical concept representing a spin- 
less particle. Here, particle concept, manifesting physical aspect of the 
Ricci scalar, is termed as riccion. It is a scalar particle with (mass)^ 
inversely proportinal to the gravitational constant. Many interesting 
consequences of dual role of the Ricci scalar are discussed here. It 
causes inflationary scenario in the early universe without taking an 
another scalar like "infiaton" also. It is found that a riccion behaves 
like an instanton also. This feature inspires "primordial inflation" . It 
is interesting to see that a riccion, obtained from higher- dimensional 

space-time, decouples into fermion and anti-fermion pair if parity is 
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violated. One-loop renormalization of riccion indicates fractal geom- 
etry at high energy. Homogeneous and inhomogeneous models of the 
early universe are derived using dual roles of the Ricci scalar. Pro- 
duction of spinless and spin- 1/2 particles, due to riccion coupling, is 
discussed here. Contribution of riccion to the cosmic dark energy is 
obtained here through one-loop renormalization and it is shown that 
dark enargy decays to dark matter during expansion of the universe. It 
inspires a new cosmological scenario consistent with observational evi- 
dences. Key Words: Ricci scalar, higher-drivative gravity, instanton, 
one-loop renormalization, homogeneous and inhomogeneous cosmolog- 
ical models, dark energy, dark matter, nucleosynthesis and structure 
formation. 



1. Introduction 

Einstein's Equivalence Principle (EEP)and General Covariance (GC) 
are two basic principles of the General Relativity (GR). According to 
EEP, in the sufficiently small region, laws of motion can be described 
in the same form as these are in non- accelerated systems free from 
gravitational field. It means that, in the sufficiently small region, laws 
of the Nature can be understood through Special Relativity(SR) in the 
framework of Minkowskian space-time. This concept is analogous to 
the concept of manifolds which are locally cuclidcan geometrical struc- 
tures. This analogy provided an insight to understand gravity through 
the geometry of the space-time. As a result, application of the riemann- 
ian geometry appeared in the arena of the physical phenomenon like 
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gravitation. Thus validity of EEP implies manifestation of gravitation 
through metrizable and differentiable manifolds of the curved space- 
time. So, using geometry of the curved space-time, Einstein could 
generahze his theory of SR to GR as he found that even a free particle 
gets accelerated there. It means that the geometry induces a physical 
force which can not be eliminated by any means, as it is due to intrinsic 
property of the space-time. Moreover, this force vanishes locally, but it 
is realized at large scales. Gravity also has similar properties. So, this 
geometrical force was identified with gravity. On this basis, Einstein 
claimed that geometry manifests gravitation. As a result, contrary to 
Newtonian theory , Einstein theory of gravity emerged as a non-linear 
theory different from other forces in Nature. 

The action for the Einstein's theory of gravitation is constructed 
using EEP and GC , where lagrangian density is R/lSnG with it! be- 
ing the Ricci scalar and G the gravitational constant. R depends on 
components of the metric tensor as well as its first and second order 
derivatives with respect to space-time coordinates. Thus the metric 
tensor is the key ingredient of Einstein's theory of gravity. 

At the classical level, this theory is in good accord with the exper- 
imental results. But, at quantum level, it is not problem - free. The 
main problem is non - renormalzability of this theory. So the need to 
construct a theory, different from GR, was realized. During the last 
three decades, different attempts were made to construct renormaliz- 
able quantum gravity. Among these attempts, higher-derivative gravity 
is an important candidate. Its action is obtained by adding higher-order 
terms of the curvature tensor like i?^, R^^W^", R^yp^jR^^^'^'^ , Di?, R^ etc. 
to the Einstein-Hilbert lagrangian R/IGttG ( where R^i, and R/xupa- are 
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components of Ricci and Riemann curvature tensor respectively). Here 
□ is the d'Alembertian operator in curved space-time defined as 



witfi g,j,„ being components of tlie metric tensor defined tlirougli tlie 
distance function, in curved space-time as dS^ — gfj,i,dx'^dx''{ii, v — 
0, 1, 2, 3). Here g is tfie determinant of gij_i,. Higher-derivative terms are 
also relevant in string inspired theories e.g. brane-gravity and Gauss- 
Bonnet theory [1]. 

The modified action contains higher-order derivatives of the metric 
tensor which is consistent with EEP and GC . At the quantum level, 
g^ji^v are supposed to be spin-2 gravitational field, called gravitons. A 
multiplicatively renormalizable and asymptotically free theory of grav- 
ity was obtained from the modified action [2, 3, 4], but it is not unitary. 
On this ground , the theory of higher- derivative gravity is rejected. 

Now the question arises whether it is the end of the story. Definitely, it 
is not so. Investigations, from time to time, have shown that there is an 
another aspect of the higher-derivative gravity which is very interesting 
and problem- free. 

The main ingredient in the action of higher-derivative gravity is the 
Ricci scalar R. It is been found that R behaves like a spinless physical 
field also in addition to its usual role as a geometrical field. In what fol- 
lows, we review important works of physicists showing physical aspect 
of the Ricci scalar manifested by higher-derivative gravity. 

In 1967,B.S.DeWitt [5] and in 1968, Zel'dovich & Novikov [6] con- 
jectured that inclusion of higher-derivative terms, in the gravitational 
action, could avoid singularity in cosmological solutions. This idea 
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was followed by Zel'dovich and Pitaevski and , in 1971, they dis- 
cussed possibility of avoiding singularity in the early universe model 
using higher-derivative gravity [7]. In 1979, Gurovich and Starobinski 
pointed out that if an arbitrary term f{R), subject to the condition 
limR-^of{R) /-R = 0, is added to the Einstein- Hilbert lagrangian, ghost 
terms do not appear [8]. Earlier to this work, Stelle [2] had also men- 
tioned that, contrary to the square of Weyl tensor 

C^iupaC^"'''^ — Rn^paR'^"'^'^ — 2Rfj,„R'^'' + (l/3)i?^, 

introduction of R'^ or f{R) does not cause the ghost problem . In 1980, 
Starobinsky proposed independently that if the sign of the R'^ term 
in the higher-derivative gravitational action is chosen properly, one 
could obtain only one scalar particle with positive energy and positive 
squared mass. He called it as "scalaron" [9, 10]. In 1984, using the 
lagrangian density, L — R-\-eR^, B.Whitt demonstrated that Einstein's 
gravity could be described by modified metric tensor components g^i, — 
(1 -|- 2eR)g^i, and the scalar curvature R with the minimal coupling to 
gravity, given by equations [11], 

R/xi^ — -^gtiuR — 87rG'7)j,y(i?), 

where 



fUR) 



6r^ 



87rG'(l + 2eRY 
In this case, action for the scalar field R is 



d^Rd^R - ~g^,{{l/2)d'^Rd,R + 



S{R) = / d^xQe^^/^{l + 2eR)- 



r 

d^Rd^R^— . 
L 6e J 
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In these equations, the scalar R behaves hke a scalar physical field 
with ^ as squared mass. Here e is the coupling constant with squared 
length dimension. The same year, Hawking and Luttrell also found the 
similar result [12, 13]. Two years later, using the lagrangian density 
like Whitt in the conformally invariant flat space-time, M.B.Mijic et 
al [14] had shown that the scalar curvature R behaved like a damped 
harmonic oscillator with the restoring force given by ^. Further, they 
have discussed that R also drives inflation when e is large. In the same 
model, it is demonstrated that, after inflationary phase, the universe 
passes through the period when R oscillates rapidly causing particle - 
production. These particles reheat the universe on annihilation [15]. 

In higher - derivative gravitational models, discussed above, either 
the gravitational constant G is taken equal to 1 or lagrangian density 
is taken as 



^ -{R + eR^). 



SttG 

As a result, (mass)^ for R does not contain the term IdnG. But, in 
a gravitational theory, role of G is crucial. So, in [16], lagrangian for 
higher-derivative gravity was taken as 



-|- higher — derivative terms 



IQttG 

giving proper role to G and the physical aspect of R. It is found that 
trace of fleld equations from this action shows that the Ricci scalar 
behaves like a physical field with squared mass dependent on G. Here, 
it is important to mention that the physical aspect of R is manifested 
only when higher - derivative terms have significant presence in the 
gravitational action, compared to j^q- It is possible at sufficiently 
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high energy. It means that physical aspect of the Ricci scalar can be 
realized at high energy. Thus , it is reasonable to infer that, in the 
high energy regime, the Ricci scalar shows dual role (i) like a spinless 
physical field as well as (ii) like a geometrical field. 

In a field theory, field is regarded as a mathematical concept as- 
sociated with a physical particle. So, if the Ricci scalar behaves as 
a physical field at high energy, there should be a particle manifesting 
this behaviour. Here onwards, this particle is called as riccion, which is 
different from scalar mode of graviton. The difference between riccion 
and graviton is highlighted in Appendix A. 

Natural units, defined as /t^ = = c = 1 ( where is Boltzman's 
constant, h is Planck's constant divided by 2ii and c is the speed of 
light) are used here with GeV is used as a fundamental unit such that 
IGeV = 1.16 X lO^^K = 1.78 x IQ-^^^m , IGeV"^ = 1.97 x IQ-^^cm = 
6.58 X 10-25 sec. 

2. A toy model to demonstrate the Dual nature 



As mentioned above, dual nature of the Ricci scalar is manifested by 
higher-derivative gravitational action. In this section, for simplicity, 
we take the following action for - gravity to demonstrate the same. 
Here, R"^ - gravitational action is taken as [16, 17, 18, 19, 20, 21, 22, 



where a and (3 are dimensionless coupling constants. The Newtonian 
gravitational constant G — Mp^ {Mp is the Planck mass). 



of the Ricci scaleir 



23, 24, 25] 




(2.1) 
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The invariancc of Sg , given by eq.(2.1) , under transformations 
9ni^ g^u + Sg/^i, yields the gravitational field equations 

—g^uOR + 2R1Rp, - ^^g^.R'^Rp,) + P{2R,^,, - 2g^,nR 

~g^uR' + 2RR^,)^0, (2.2) 

where {; } denotes covariant derivative in curved space - time and □ is 
defined by eq.(l.l). 

Taking trace of equations (2.2), one obtains 



OR + m'^R = 0, (2.3a) 

where 

- [327r(a + 3^)^]-^ (2.36) 

This equation shows that R can be treated as a mathematical concept 
representing a spinless particle with squared mass , given by eq.(2.3b). 
To understand the physical meaning of the equation (2.3), it is useful 
to take an analogy of the Klein - Gordon equation for a scalar field 

□0 + mJ(/) = O. (2.4) 

Equation (2.3) seems to be analogous to the equation (2.4), but there 
are two problems : 

(i)Mass dimension of is 1 as eq.(2.4) is derived from the action 

^^ = IJ d'xy^[g^''d,(l>d,(l> - my]. (2.5) 

But mass dimension of it! is 2, being combination of square of first 
derivatives and second derivative of components of the metric tensor. 
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Qiiu {g^iu are dimensionless as dS^ = g^i,dx^dx" .) with respect to space 
- time co-ordinates. 

(ii) □ and both depend on g^j^i,, but (j) does not depend on it. 

In what follows, these problems are resolved. First problem is solved 
multiplying the equation (2.3a) by 77 ( where , in natural units, \r]\ — 1 
with {massy^) and recognizing rjR as R. As a result, the equation 
(2.3a) is rewritten as 

D^ + m^^^O (2.6) 

It is easy to find that mass dimension of i? is 1 as mass dimension of 
i? is 2 (mentioned above). 

The second problem is resolved through the following discussion. In 
the small vicinity of a space - time point, g^ui^) can be expanded as 

9fiu{x) = 6/^^(0) + ^R^au)3{0)x'^x'^ - ^d^R^aup{0)x'^xf^X^ + 

(2.7) 

where g'^^(O) = 77^^. 

Using this expansion, one can easily see that the Ricci scalar R is 
non-zero in a locally inertial co-ordinate system, unless the entire space- 
time manifold is Minkowskian. Also, in the locally inertial co-ordinate 
system, □ operator (defined by eq.(2.2b)) looks like 

□ = ri^" , 

^x^^^x'^ 

showing that, in a locally inertial co-ordinate system, □ is a similar 
operator for R as it is for other scalar field 0. According to the principle 
of equivalence and principle of covariance, this characteristic feature of 
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□ and R (these being scalars) will be maintained at global scales also 
[25]. 

On the basis of these analyses, it is obtained that the Ricci scalar 
not only behaves as a geometrical field, but also as a spinless matter 
field analogous to (f). 

To exploit the matter aspect of the Ricci scalar, R is treated as a basic 
field. So there should be an action Sp^ yielding eq.(2.6)as a result of 
invariance of 6"^ under transformation i? — > R-\-5Rva the same manner 
as Sfj, yields eq.(2.4)using invariance under transformation — > + 50. 

If such an action 5"^ exists, 

SS^^- J d^x^ SRp + m^]R, (2.8) 
implying eq.(2.7) subject to the condition SSp^ = 0. Prom eq.(2.8). 



-I 



d'^Xy 



-9 

-g 6 



g^^d^.Rd^idR) - rri^RdR^ 



^"d.Rd^R - m'R 



(2.9) 



In a locally inertial co-ordinate system. 



J d'^xy^ 5| ^^d^Rd^R - m^R'j } 
= J d^X 5| ^^df^Rd^R - m^lfj } 



,2 ^ 



} 



as i?,being a scalar, is equivalent at local as well as global scales. Here 
X*(i = 0, 1, 2, 3) are local and a;*(i = 0, 1, 2, 3) are global coordinates. 
The action 
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being covariant, is invariant under transformations X* — > (i.e. from 
local to global coordinates), so 



Thus, eq.(2.9) reduces to 



^d^Rd^R-m^R^ } 



-gf^'d^Rd^R - m^fej } 



yielding 



d^x^/^ ^g^'^df.Rd^R - m^R^ 



(2.10) 



Thus eq.(2.10) gives the required action for the spinless physical field 
R. It is different from scalaron, proosed by Starobinsky [9, 10] in two 
ways (i) mass dimension oiriccion is 1, whereas the same for a scalaron 
is 2 and (ii) (mass)^ for a riccion is inversely proportional to G (the 
gravitational constant), but a scalaron is independent of G. 

3. Riccion as an instanton 

Here gravitational action is taken as [22] 



-9 



R 



IGttG 



- aRmuuR"" - PR^ + WiR" + 9ni?')], (3.1) 



where a, (3 and A are dimensionless coupling constants. G and 7] are 
defined in the previous section. Though term containing \^R? is a 
surface term, still it is needed here. The motivation for this type of 
action has been clearified below. 
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The action (3.1) yields gravitational field equations as 



IGttG 
1 



-g^r^RP'^R,,) - (3{2R,^, - 2g^,UR - -g^^R' + 2RR^,) 
+Xr]''{12Rl, - ^g^.OR' + ISnR" - ^-g^^R^ + SR'R^,) = 0. 



(3.2) 



(3.3a) 



Trace of equations (3.2) is obtained as 

OR - m^R + Xif{a + 3(3)'^R^ = 0, 
where 

= [327r(a + 3P)G]-\ (3.36) 

Multiplying eq.(3.3) by rj and recognizing rjR as i?,eq.(3.3) is re- 
written as 



HR^ ^r(r'~^) (3.4) 

Proceeding on the line of the previous section, action for riccion is 



obtained as 



-I 



-g 



-g^'^di.mR-ViR) 



with 



(3.5a) 



^ ' 8(a + 3/5)V 327rA/ 
Lagrangian, obtained from eq.(3.5), is 



(3.56) 



L{t) = / d'^x^ -g^^d^Rd^R - V{R) 



with V{R) given by eq.(3.5b). The corresponding hamiltonian [22] is 



rl 



H{t) = E{t) = / d-'x^ -g^'d^Rd.R + V{R) 



(3.6) 
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where integration is performed over t = constant hyper surface. As the 
potential V{R) does not contain the velocity like term R, the hamil- 
tonian gives the total energy E of riccion R in this case. 

Now, one can find the condition for the finite euclidean action cor- 
responding to the lagrangian density for R, which gives the classical 
path for imaginary time called the instanton [26, 27]. 

If the early universe is spatially flat, homogeneous and isotropic, its 
geometry is given by the line-element 

dS"^ ^ dt^ - a\t)[dx^ + dy^ + dz^] (3.7) 

In the background geometry of this model, the riccion energy (given 
by eq.(3.6)) can be written as 

where v is the volume of the 3-dimensional space and V{R) is the 
potential,given by eq.(3.5b). 

Now, using the transformation t —iT, one obtains that the Ricci 
scalar R transforms to —R. As a result. 



R^-R (3.9) 

and 

^-«'wH-i(f)'+''(^)}- ("°) 

Thus, in imaginary time, R can go from the state R = , = to 

^ y/327T\{a+3f3) 

the state R = + , . with ^ = 0. So, setting ^ = in eq.(3.8), 

one obtains the classical trajectory in the R — V{R) with imaginary 
time as 
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Mp , / MpT \ 
R= J— tanh( — 3.11 

which is an instanton solution. If couphng constants A, a and /3 satisfy 
the equation 



A((5 + 3/3)-l, (3.12) 
the action for this trajectory, in the imaginary time system, can be 
calculated as 

(3.13) 

using eq.(3.10) with = 0. In what follows, it will be shown that Se 
is finite, which is essential for riccion to behave like an instanton. 

The basic equation showing the physical aspect of i?, in this case, 
is eq.(3.4). So, i?, given by eq.(3.11), should satisfy it for the sake of 
consistency. In the background geometry of the early universe,given 
by eq.(3.7), eq.(3.4) is re- written as 

~ 3d ~ A^ ~ /~2 Ml\ , 

Connecting eqs.(3.13) and (3.14), one finds that eq.(3.14) is satisfied 
by given by eq. (3. 11), only when 

> ^I:^ ^ 2.45 X 10-^^GeV-i. (3.15) 
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as tanhl2.2 = 1 and sechl2.2 = 0. For t^/\ < = 2.45 x 

10~^''GeV~^, eq.(3.15) yields ^ = 0, which satisfies eq.(3.13),only when 
T — 0. So, to be on a safer side, one can take ry/X > 2.45 x 10~^^GeV~^. 
As a result, eq.(3.11) yields 

R=^^. (3.16) 
4^2^ 

With t — —IT, R is obtained from the line-element, given by eq.(3.7) 

as 

R^-qT-+(^)\ (3.17) 
La \ a/ J 

where prime denotes differentiation with respect to r. Thus, for r\/A > 
2.45 X 10-i^GeV"\ one obtains from eqs.(3.16) and (3.17) 

(3.18) 

This equation yields the solution 

"°°'4Kiife)""(^-^»'] p-"") 

with 



24r7V2^\i/2 



Mp 

Connecting eqs.(3.13) and (3.19), one finds that 



(3.196) 



- A V Mp J V4V2^/ /T/ 



= 0.62 X tt-'/'X'/'m'J'I^ 
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Using T = it in eq.(3.19a), wc get 
a{t) = aocosh^'^ h (— 



Mp 



) 



1/2 



{t - to) 



(3.20) 



It is discussed above that eq.(3.14) is valid when t = |r| > 2.45 x 
10-^^GeV"\ so we can take to = 2.45 x IQ-^^GeV"^ = 1.61 x IQ-^^sec. 
in eq.(3.20). Thus to is well before the grand unified phase transition 
which is expected to occur around 10~^^sec.. This phenamenon is called 
primordial infiation [28, 29]. 

. Thus, it is seen riccion behaves as an instanton too and the in- 
stanton solution, so obtained, inspires inflationary scenario which starts 
much earher to GUT phase transition. 

4. Emergence of riccions and riccinos from higher - derivative 
higher - dimensional gravity 

Here, the line-element for (4 + D) -dimensional space-time is taken 
as [20, 30] 

ds^ = Qf^^dx^'dx'' - pldel - pldel pidel, (4.1) 

where i/ = 0, 1, 2, 3 and yi = piOi, y2 = ^26*2, ■ ■ • ,yD ^ PdOd{^ < 
^1, ^2, • ■ ■ < 27r) as is topologically equivalent to product of D 
copies of circles with different radii pi, P2, ■ ■ ■ ^ Pd- 

In this case, the gravitational action for higher- dimensional higher- 
derivative gravity is taken as [20] 



where = (27r)^pip2 • • • PdG, is the determinant of qmn 

{M,N ^ 0,1,2,- ■■ ,{D + 3),a ^ a X [(27r) VP2 • • • Pd]~^ and P = 




(4.2a) 



DUAL NATURE OF THE RICCI SCALAR 17 

P X [(27r)^pip2 ■ ■ • PdY^ ol and fi are dimensionless coupling constants 
in the four-dimensional theory) and 



1 d 



d 



The action(4.2) yields gravitational field equations 

-{RmN — -9MnR{A+D)) + Oi{'2R;MN — '^9Mn'^{A+D)R{4:+D) 



167rG(4+D) 2' 

~2^MnR(4+d) + '^R{A+D)Rmn) + /3(2n(4+D)i?MiV 

-^fl'AfiVn(4+D)%+D)) = 0, (4.3) 

Prom the line-element, given by eq.(4.1), i?(4+D) = i?4. So, taking trace 
of eqs.(4.3), one obtains 



UR + XR^ + m^R = 0, (4.4a) 

where i?4 = i?, ^4 = □. 



(4 46) 

167rG'[4(L> + 3)a + Dp] ^ ' 



and 



A = ^ (44c) 

167rG[4(D + 3)a + ^ ' ' 

subject to the condition [A{D-\-?,)a-\-DP] > 0. This condition is imposed 
to avoid the ghost problem as above. The operator 0(4+^) effectively 
reduces to ^4 as R is independent of y- coordinates. 

Thus , from eqs. (4. 4), equation for R is obtained multiplying it by t] 
(which is used for dimensional correction as explained in section 2, as 



UR + XRR + m^R = 0, 



(4.5) 
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as a result of spontaneous compactification. It is interesting to note 
that A = for D = 0, which corresponds to the case of 4-dimensional 
space-time geometry, discussed in sections 2 and 3. So, the term XRR 
is a manifestation of higher-dimensional geometry in four-dimensional 
space-time. 

In what follows, two cases are discussed (i) when A = 1/4 (ii) when 
A 7^ 1/4. The first case is possible only when f3 i.e. in the presence 
of □(4+D)-R(4+D) in Sg, given by eq.(4.1). In this case, eq.(4.5) looks 
like 

(n + ^R + m'^^R = 0, (4.6) 

with given in eq.(4.4a). 

One can express the scalar field R as 'ipip ( where -0 is a Dirac spinor, 
ip = -j^t^o ^[^i^ as hermitian conjugate and 7° as a Dirac matrix in 
flat space-time). Physically, it means that the riccion can be envisaged 
as a composite of spin-1/2 fermion and anti-fermion [31]. For this, it 
is necessary for eq.(4.6) to decouple into Dirac equations for and ■0, 
which is possible under a condition derived below. 

Equation (4.6) can be written in the matrix form as 



AB^Z (4.7a) 

with 

A^(a + ^R + m^y, (4.76) 

B = RI (4.7c) 

and Z is the zero matrix which is a 4x4 matrix with all elements 0. / 
is the 4x4 unit matrix. 
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The operator (n + ^R^I can be written as j^V/^Y^u [32, 33, 34]. 
So, the matrix operator A looks hke 

A - (t'^V^t'V. + m'/) (4.7d) 

Here are Dirac matrices in curved space-time related to 7", Dirac 
matrices in flat space-time as 

r = • (4.8) 
7^^ and 7" satisfy anticommutation rules 



{7^7'^} = 2^'^'^ (4.9a) 

and 

{7", 7'} = 277'^^ (4.96) 
where 77"* (a, 6 = 0, 1, 2, 3) are components of Minkowskian metric ten- 
sor, e^, in eq.(4.8) are tetrad components defined as 

g,^^ = e^e^77«^ (4.10) 

The operator (covariant derivative in curved space - time), given 
above, is defined as 



= 9^ - r^, (4.11a) 

where 

= -\{d,e^a + K^el)g^,elfr (4.116) 
with being Riemann-Christoffel's symbols. 
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Using R = ipip in eq.(4.7a) and algebra of Dirac matrices, one obtains 



V' (7^ - im) (7'' + im) V' = (4. 12a) 

provided that 

V^V.V^t'^T^ = 0. (4.12&) 
This condition can also be expressed as 

"^^NMrYi^R = -V^V^V^rt^t^l, (4.13) 

where ipL,R = [(1 ± 7^)72]^^ and 'iI)l,r = il^iil =F 7^)/2] with 7^ = 
^7°7^7^7^- The equation (4.12a) imphes the Dirac equation for ^0 as 

(7''V^ + im)^ = (4.14a) 
and its complex conjugate 

i^{YV^-im) = Q. (4.146) 

This shows that eq.(4.7a) imphes eqs.(4.14) under the condition pro- 
vided by eq. (4. 13), which shows the breaking of left- right symmetry. 
Physically, it means that a riccion (with A = 1/4) can decouple into a 
fermion and anti-fermion pair if left - right symmetry breaks [20] . In 
other words, riccion can decompose as 

R^il) + ip, (4.15a) 

provided that parity is violated. Eqs.(4.14) indicate that spin-1/2 
particles, so obtained, are electrically neutral showing that ip and ip 
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will be invariant under charge conjugation ( C-symmetry). So, CPT- 
invariance implies that the decoupling of riccion to pair of riccinos and 
anti-riccinos (fermionic and anti-fermionic counterpart of a riccion) will 
not respect T-symmetry (time- reversal) [35] . It means that riccinos and 
anti-riccinos can not combine to form riccions i.e. it is not possible to 
have back reaction 

ij + tp^R, (4.156) 
Here CP violation is taking place , so it may explain excess of baryons 
over anti-baryons if riccions decay to observed baryons. 

If A 7^ 1/4, the operator (D + Ai?) can not be expressed as 7'^V^7'^Vi,, 
so eq.(4.6) can not be written as (4.12a). So, in this case, decoupling 
of a riccion to a riccino and an anti-riccino is not possible. 

According to the above theory, it is tempting to speculate that our 
physical universe might have emerged from higher-dimensional geom- 
etry in the extreme past through decay riccions to riccions and anti- 
riccinos, when parity is violated. 

5. One - loop renormalization of riccion obtained 
from the higher- dimensional space-time 

In what follows, lagrangian density of the riccion with a quartic self 
interaction potential is obtained from higher-derivative gravitational 
action in (4-|-D)-dimensional space-time. 

The topology of the space-time is taken as (M^ as a 4- 

dimensional space-time and as D-dimensional torus as an extra- 
dimensional compact manifold) with the distance function defined as 

dS^ = g^^dx^dx'' - pldel - pldOj pldOl, (5.1) 
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where gfj,u{fi, ly = 0,1, 2, 3) are components of the metric tensor, pi{i = 
1,2,- ■■ ,D) are radii of circle components of T^[D ^ 2) and < 
9i, ^2, • ■ ■ < 2n. As usual, the space-time manifold is taken to be 
C°°-connected, Hausdorff and paracompact without boundary [22, 23, 
24, 25]. 

Here action for the higher-derivative gravity is taken as [24] 

/r R 
d}xd^y^-g(^i+D) +«(4+g)i^(4+D)+7(4+D)(i^?4+g) 

— [D - 2) (^+^)-^(^+4)^J' [p.Za) 

where G(^i+D) = GVd, «(4+i:>) = aVj^^, 7(4+d) = 3!(j_2) ^g^ 

yz? = (27r)Vp2---PD (^7^2), i?B = 

and g(4+D) is the determinant of the metric tensor 5'yv/iv(^; iV = 0, 1, 2, • • • , 
(4 -|- D). Here G(4+d) is the (4-|-D)-dim. gravitational constant and 
Q!(4+D) as well as 7(4+d) are coupling constants, o; is a dimensionless 
coupling constant, i?(4+D) is the (4-|-D)-dim. Ricci scalar and G is the 
4-dim. gravitational constant. 

From the action (5.2a), gravitational field equations are obtained as 

(l67rG(4+£)))~'^(-RMAr — -j^gMNR{A+D)) + Q;(4+D)-H"mjv + 7(4+d)-H"mjv 

(5.26) 

where 

-^MiV ~ '^R;MN — '^9MN^{4+D)R{4+D) — l^^^^^'i^+D) + '^R{4:+D)RmN , 

(5.2c) 

and 

f7-{2) _qp2 q„ n r2 6(D4_3) 1 n ^2 

^MN — '^^■,MN~'^9MN^{4+D)^(4+D)~ -^j^—^\^2^MN^{4+D)^(4+D) 
+'^'^{4+D)R{D+4)RmN + R^Mn} ~ T^QMnRIa+d) + '^R{4+D)^MN (5. 2d) 
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with semi-colon (;) denoting curved space covariant derivative and 
n 1 9 f ^- MN d \ 



Trace of these field equations is obtained as [24] 
- "^^^^^^ ~ q;(4+d)[2(-D + 3)n(4+£,)i?(4+D) + 2-^-^(4+0)] 

-1(^+D)[+\{D - 2)i?f,+^)] = 0. (5.2e) 

Using the method of earlier sections, the equation for riccion is ob- 
tained as 



[n + -^R + m^ + ^R^]R^O (5.3a) 
2 3! 



with 



D 



A = 



2(D + 3) 

647rG'a(L' + 3) 
1 



4(D + 3)q;' 

(5.36, c, d) 

where o; > to avoid the ghost problem [9] . As done in earher sections, 
the corresponding action for R is obtained as 

S^ = ^ d'x^g^-g^^d,Rd^R-{^^^iR^-r\^^ (5.4) 

5 (a). One-loop quantum correction and renormalization 

The 5"^ with the lagrangian density, given by eq.(5.4), can be ex- 
panded around the classical minimum Rq in powers of quantum fluc- 
tuation Rq = R — Rq as 



24 

where 



sf 
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/ 



d x\/—g 



A ^ 



S 
and 



f = jd^xV^R,[a + ^R + m' + ^Rl]R, 



S 



(1) 

R 







as usual, because this term contains the classical equation. 

The effective action of the theory is expanded in powers of h (with 

^ = 1) as 

r(^) ^s^ + r(^) + r' 

with one-loop correction given as 

r(i) ^ 



where 



D 



R 



5R^ 



-lnDet{D/n^), 



R=Ro 2! 



(5.5a) 



(5.56) 



and r' is a term for higher-loop quantum corrections. In eq.(5.5), /i is 
a mass parameter to keep T'^^^ dimensionless. 

To evaluate T'^^\ the operator regularization method [36] is used. 
Up to adiabatic order 4 (potentially divergent terms are expected up 
to this order only in a 4-dim. theory), one-loop correction is obtained 
as 



r(^) = (167r 



ds 



-i- 

6V5 



(s-1) 

—R'"'Ra^ + -(--n R^ \ 
180 " 2^6 J J 



s-2)(s- 1) 



s=0 



(5.6a) 
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where 



M^ = m' + {\/2)Rl. 



(5.66) 



Here it is important to note that matter as well as geometrical both 
aspects of the Ricci scalar are used in eq.(5.6). The matter aspect is 
manifested by R and the geometrical aspect by i?, Ricci tensor compo- 
nents R^j^i, and curvature tensor components R^ua/s as it is mentioned 
above also. 

After some manipulations, the lagrangian density in F^^^ is obtained 



Now the renormalized form of lagrangian density can be written as 



as 




(5.7) 



ren 



^g'^'^d.Rod^Ro - ^Rl - ^rn'Rl - + A 
+eQR+ -eiR^ + e2R^^ R^u + ^^R^^'^^ Rnuap 

+e4ni? + rw + Lct 



(5.8a) 



with bare coupling constants Aj = (m^, A, A, ^, eo, ei, 62, £3, 64), F*^^) given 
by eq.(5.7) and L^t given as 



^ct = -^^diRRl-^^drn'Rl-^-^Rt + SK + Se^R + ^^de^R^ 
^Se^R^'^R^, + Se^R^^^^R^^^p + Se^UR. 
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(5.86) 



In cq.(5.8b), 6Xi = {6m^,SX,SA,S^,Seo,Sei, 662,663, 664) are counter- 
terms, which are calculated using the following renormalization condi- 
tions [37] 



A 
A 



= L 



64 



dRt 



Ro=R^o)i,R=0 



dRl 



Ro=0,R=0 



Q3 

—77 —Lj. 

dRdRl 
dR ''''' 



d 



Ro=0,R=0 
_Ro=0,i?=i?5 



d{R^^^R^,) 
d 



Rq=Q,R=Rq 
Lrf 



Rq=o,r=R8 



(5.9a, 6, c, d,e,f,g,h,i) 



As R — rjR, so when R — 0, i?(o)o = R{q)i — R(q)2 — and R5 
i?6 = i?7 = i?8 = when ^0 = 0. 



DUAL NATURE OF THE RICCI SCALAR 27 

Eqs.(5.8) and (5.9) yield counter-terms as 

IGn'SC = -3A(e-^)/n(mV/x2) 

16 = m2(^C- i)/n(mV/x2) 

167r^5e3 = Y^/n(mV/^^) 

16 A4 = ^(^-e)/n(mV/.2) 



(5.10a, 6, c, d,e,f,g,h,i) 



(b) Renormalization group equations and their solutions 



The effective renormalized lagrangian can be improved further by 
solving renormalization group equations for coupling constants Xi^R) (suffix R 
stands for renormalization, which is dropped onwards). For this pur- 
pose one-loop /9-functions, defined by the equation [32, 37, 38] 

Px,^pi-^{K + 5Xi)\^ (5.11) 
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with counter-terms 6Xi from eqs.(5.10), are obtained as 



Pa 

Pe, 
Pes 



4 

m 



167r2 
3A2 



87r2 



87r2 

3A(^-i) 



87r2 

Jill) 

87r2 
87r2 

1 



14407r2 
1 



14407r2 



(5.12a, 6, c, d,e,f,g,h,i) 



using the fact that = for bare couphng constants Aj 

The renormahzation group equations are given as 



(5.13) 



where t — ^ln{rn^^/ ^'^) with being a mass parameter and rUc being 
a reference mass scale such that /i > rric- Using /3-f unctions for dif- 
ferent couphng constants given by eqs.(5.12), solutions of differential 
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equations (5.13) are derived as 

,4 
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A 
A 



■rUn 



— ^^0 



3Aot 



3Aot^l/3 
8^ 



m 



Ao 

ml 

1 
6 



1 - 

1 - 

'^0 



87r2 
3Aot 



-1/3 



3Aot 
8^, 



-1 



64 



^00 + 

eio + 

£20 + 
^30 — 

£40 — 



Ao L V Stt^J \ 



t 

14407r2 
14407r2 

14407r2 







18An 



3Ant^ 



(5.14a, 6, c, d,e,f,g,h,i) 



where Ajo = Aj(i = 0) and t — at /i — rric according to defini- 
tion of t given above. 

These results show that as ^ oo(t — > — oo) A ^ and 
and ^ i. Thus it follows from these expressions that in the limit 
/X — > cxD, theory is asymptotically free. 

Further it is assumed that D = at energy mass scale /i — rric- Now 
recaUing the definition of ^ from eq.(5.3b), one obtains ^o = and 



D 



1 /-I SApt 
^ \^ 87r2 



2 + 



(5.15) 
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It is interesting to sec from this result that D increases with increas- 
ing energy mass scale ji and it is equal to 1 when = 1.8 x l^^^ruc 
provided that Aq = 1. Moreover, D is not necessarily an integer but a 
real number. This result suggests that dimension of the space-time will 
be equal to 4 at = mc and will increase continuously with increasing 
energy mass scale. Non-integer vahics of the dimension indicate that 
the space-time above the energy mass scale = iric should be frac- 
tal [39, 40], as fractal dimensions need not be integers like topological 
dimensions. 

Eqs.(5.3d) and (5.14b) yield 



a — 



1 r 
12 



t 



47r2 



(5.16) 



and QJo = if-^o = 1- 

Recalling the definition of from eq.(2.8c) and using eqs.(5.15)- 
(5.16) with Ao = 1, one obtains 



G = G 



N- 



6 



47r^ 



1 _ 

^ 87r2 



2/3' 



(5.17) 



where Go = Gn{Gn is the Newtonian gravitational constant). Now 
from eq.(5.3c) 



1 



rrin 



Ml 



(5.18) 



° 87rGiv 87r 

as = at t = or/i = rUc and Gn = Mp'^{Mp is the Planck's 
constant). Thus eqs.(5.14a) and (5.18) imply 

A-Ao 



2(87r)^ 



87r2. 



(5.19) 



The equation (5.19) implies that vacuum energy density A increases 
with increasing energy mass scale /i ( or decreasing t ). Taking Aq = 1, 
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it is obtained that (A — Aq) — > oo when /i oo{t — > — oo). In the 
following table (A — Aq) is exhibited for different values of /j,. 

Table no.l 

(A — Aq), G/Gn and dimension of the space — time (4 + D) are 
tabulated below against /i/mc with mc = 1.76 x lO^^GeV taking Aq = 



ji/mc 


(A - Ao) 


G/Gm 


(4 + L>) 




in GeV^ 






1 





1 


4 


1 + 10-6 


9.58 X 10^^ 


1 


4 + 3.6 X 10-^ 


1 + 2 X 10- 


-6 2.04 X 10*^5 


1 


4 + 7.7X 10-8 


1 + 3 X 10- 


-6 2.99 X 10^5 


1 


4 + 1.1 X 10-^ 


1 + 4 X 10" 


-6 4.07 X 10^5 


1 


4 + 1.5 X 10-^ 


1 + 10-5 


1.01 X 10*56 


0.9999999 


4 + 3.8 X 10-^ 


1 + 10-^ 


1.01 X 10^^ 


0.9999998 


4 + 3.8 X 10-6 


1 + 10-3 


1.004 X 10^8 


0.9999993 


4 + 3.8 X 10-5 


1 + 10-2 


9.995 X 10^8 


0.9999389 


4 + 3.78 X 10-4 


1.1 


9.56 X 10^9 


0.9993944 


4 + 3.6 X 10-3 


2 


6.903 X 10™ 


0.9995533 


4 + 2.6 X 10-2 


10 


2.23 X 10^^ 


0.9846903 


4 + 8.2 X 10-2 


10^ 


1.02 X 10^2 


0.9199941 


4.339 


10^° 


1.998 X 10^2 


0.825 


4.587 


1020 


3.3 X 10^2 


0.705 


4.826 


1030 


4.35 X 10^2 


0.618 


4.966 


oo 


oo 





5.5 



Here, it is interesting to note that when ji comes down from infinity to 
l.OOOOOlmc , (A — Ao) decreases slowly. But when it comes down from 
l.OOOOOlmc to mc, there is a sudden drop in the value of (A — Aq) from 
9.58 X lO^^GeV^ to zero. It means that all of a sudden, a huge amount 
of energy ( with density equal to 9.58 x lO^^GeV^ ) is released Bi ji — 
l.OOOOOlmc, which is sufficient to heat the universe up to a temperature 
1.76 X 10^6QgV. It corresponds to the energy scale 1.76 x lO^^Qgy 
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in natural units ( as Boltzmann constant /t^ = 1 in these units ). 
Sudden release of energy indicates a phase transition at 1.76 x lO^^GeV. 
According to the standard model of grand unified theories, strong and 
electro-weak interactions unite at lO^^GeV. But gravity maintains its 
identity different from these interactions up to this energy scale also. 
So unification of these interactions are expected above lO^'^^GcV [31]. 
Thus symmetry breaking due to phase transition at 1.76 x lO^^GeV 
( discussed above ) is expected to be the energy mass scale where 
gravity decouples from strong-electroweak interactions. In other words, 
these results suggest that unification of gravity with strong-electroweak 
interactions should take place at 1.76 x lO^^GeV. It implies that 

rric ~ 1.76 x lO^^GeV. (5.20) 

As discussed in the preceding section, another interesting result is 
given by eq.(5.15) which suggests that above the energy mass scale /i — 
l.OOOOOlmc ~ 1.76x 10^*^GeV, space-time should be fractal . Moreover, 
according to this result, dimension of the space-time increases from 4 
to 5.5 with the increasing mass scale. 

The equation (5.17) imphes that the gravitational constant G de- 
creases with increasing // with G = Gat at // = 1.76 x lO^^GeV. 

Eqs.(5.14) also show that a, eo, ei, €2, €3 and £4 increase with increas- 
ing n showing that higher- derivative terms grow stronger with increas- 
ing energy mass scale. 

6. Cosmological models inspired by dual role of the Ricci Scalar 

and peirticle production 

6.1. Homogeneous as well as inhomogeneous 
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models of the early universe 

Here, the gravitational action with thermal correction is taken as 
[21, 23] 

Sg = J - + ^XT\a + 3P)R + aR^uuR^''+ 

pR^ - 2X7]\d + 3(3) {R^ + 9aR^)] , (6.1.1) 

where a, $ and A are dimensionless coupling constants. Here T is the 
temperature as well as rj is another constant with length dimension 
used for dimensional correction [16, 17, 19, 20, 21, 22, 23, 24, 25], as 
given in earlier sections. 

Sg ( given by eq.(6.1.1)) yields gravitational field equations 



^hT\a + 3/3) - ^ {R,, - '-g,,R) + a(2i?^^,^ - UR^, 



1 

T 



^g^^OR + 2R1Rp, - ^g^^R'^'^R,,) + P{2R.^, - 2g^,nR 
+ 2RR^,) - 2\rf{a + 3;3)(12i?f^, - ^g^^UR^ 



+18ni?' - ^^g^.R' + 3R^R^,) = 0, (6.1.2) 
Trace of these field equations is obtained as 

UR = {m" - ^XT^)R^ - Xrj^R^+, (6.1.3) 

where 

= [327r(a + 3(3)G]-^ (6.1.36) 
Using the method ,used in section 4, eq.(2.3) is re- written as 



BR = (m^ - ^XT^)R - XrfR^ (6.1.4) 

Like earlier sections, here also this equation is multiplied by 77 and 
■qR is recognized as ^. As a result,the equation (6.1.4) looks like 
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1 



OR = {mf - -XT')R - Xtj'R'^). (6.1.5) 
The action for R, yielding the equation (6.1.5), is obtained here, 



following the method of earher sections, as 



S 



R 



(6.1.6a) 



where 



V^{R) = --m^R^ + —T^) + -XR^ + -XT^R'' - —T\ (6.1.66) 
^ ^ 2 ^ 12 ^ 4 8 90 ^ ' 

Spontaneous symmetry breaking can be discussed with the Higgs - 
like potential V^{R) given by eq.(6.1.6b). Vacuum states are obtained 
from [16, 23, 31] 



OR 



0, 



or 



-m^R + XR^ + -XT^R = 0, 
4 



yielding turning points of V'^{R) as 



and 



where 



R^O 



2m 

7t 



(6.1.7) 



(6.1.8a) 



(6.1.86) 



(6.1.8c) 



These states can be written in terms of vacuum expectation value 
of ^ as < ^ >= and < ^ >= (1/2) ^(r^ - T^). It shows that as 
long as T > Tc , the field R remains confined to the state < R >— 
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0. But when temperature goes down such that T Tc , R tunnels 
through the temperature barrier T = Tc and settles in the state < 
R >= (1/2)7; when T < T^. The state < ^ >= is called the false 
vacuum and < R >— (1/2)7; corresponds to the true vacuum state . 
Moreover , symmetry is intact in the state < i? 0, but it is broken 
spontaneously at < i? >= ±(l/2)Tc. As a result, energy is released as 
latent heat in the form of radiation with density given as 

V^(< i? >= 0) - V^{< R >= ^Te) = ^. (6.1.8d) 

(a) Homogeneous models of the universe 

In the homogeneous case, the line-element of the cosmological model 
of the early universe can be taken as [16] 



where a{t) is the scale factor and t is the cosmic time. The 4-components 
of a time-like vector u'^ obey the normalization condition 



uf'uu = +1. 



(6.1.10) 



For the model, given by eq.(6.1.9), it! is obtained as 

R=6r) 



'd 




- + 




La 





(6.1.11) 

where dot(.) stands for derivative with respect to the cosmic time t. 
In the symmetric state < R >— 0, when T > the equation (6.1.11) 
reduces to 



a (a 
— + 
a 



a / 



0, 



(6.1.12) 
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where prime (/) means derivation with respect to {t/tp), with tp being 
Planck's time. 

Solution of the equation (6.1.12) is obtained as 



a^(t)^a^o + ^ (6.1.13) 



1^1 
tp 

with a{t = 0) = flfl. It shows that the model expands as (og + |^|/^p)^^^ 
in the symmetric state above the critical temperature Tc. Now, it is 
important to decide whether qq — at t — 0. Answer to this question 
is obtained as follows. 

Components of the energy-momentum tensor for riccion , obtained 
from the action Sj^, given by eq.(6.1.6), are defined as 



2 5S^ 



yielding [16, 23] 



2r)R.,DR 



+2r}m^RR^^ - 2r]XR^R^^ - ^rjXT^RR^^ + Ar^m^R.^^ - irjm^g^^nR 
-4:riXRl, + Ar]\g^,UR - XT^R^, + XT^g^^nR. 

(6.1.14) 

Using eqs.(6.1.10) and (6.1.14), one obtains that 



Ti.^u^u'' < (6.1.15) 
in the state < R >= 0, implying violation of the energy condition at 
and above the critical temperature. This result shows that the model 
will bounce at t — o without an encounter with singularity. As a reult, 
it is obtained that 
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Oo 7^ 0. (6.1.16) 

Due to expansion of the model, given by eq.(6.1.13), temperature 
falls down. As discussed above, when it is sufficiently below the critical 
temperature such that <^ T^^, the model acquires the state 

< ^ >= (6.1.17) 

Connecting eqs.(6.1.11) and (6.1.17), a differential equation is ob- 
tained as 



a /d\2 Tr 



a \aJ 



1277 



yielding the solution 



(6.1.18) 



a^{t) — \l ——sinh 

^-^ c 



24 




al [sinhD] ^sinh (t — tc) \ I — + D 

67] 




a^exp 



(6.1.19) 

where Oc and tc are constants such that Qc — a{t — tc) with tc, being 
the time when this type of expansion of the model begins as well as 

D — sinh~^ 



TjTc 



The second case is 



a \aJ 



24 



T 

c 

12^ 



(6.1.20) 
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yields the solution 




(6.1.21) 



where Oc and tc are constants such that Gc = a{t = tc) with tc, being 
the time when expansion of the model follows eq.(6.1.21). 

But a{t), given by eqs.(6.1.19), shows that the cosmological model 
will expand more rapidly in the state R —< R >= (l/2)Tc compared 
to the state R =< R >= 0. The expansion , in the state R =< R >= 
— (l/2)Tc, is given by eq.(6.1.21)., which shows that — a^g < < afg 
leading to imaginary a{t) also, which is unphysical. So the state < 
R >— — (l/2)Tc will not be considered any more. 

(b) Inhomogeneous models of the universe 

First type of inhomogeneous models 

Though large-scale observations, in the present epoch, support ho- 
mogeneous and isotropic cosmological models, still it is believed that a 
realistic model of the universe should not be exactly homogeneous, at 
least, at small scales. This idea motivates to investigate cosmological 
models, being spatially inhomogeneous and anisotropic at small scales 
but asymptotically homogeneous and isotropic. 

Spatially homogeneous and anisotropic cosmological models such as 
Kasner model, Lukash and Starobinsky, Bianchi-typel and type IX 
models [34] have been proposed from time to time in the past. But more 
complicated spatially inhomogeneous and anisotropic models were in- 
vestigated during the last two decades [41, 42, 43, 44, 45, 46]. Due to 
mathematical difficulty, in solving Einstein's field equations, even these 
models are inhomogeneous with respect to only one spatial coordinate. 
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Here, cosmological models, being inhomogeneous and anisotropic 
with respect to all three spatial co-ordinates in the (l+3)-dimensional 
space-time, are derived using the method employed in section 6.1(a). 

The general framework of the line-element, for cosmological models 
of the early universe, is taken as [21] 



where f{x, y, z) is the function resposible for the inhomogeneity. 

Looking at the line-element (given by eq.(6.22)), the reader may 
think that the function f{x,y,z) can be absorbed redefining coordi- 
nates {x, y, z) as new coordinates (X, F, Z) such that 



At this stage, it is important to mention that redefinition of coordi- 
nates is not valid at points, where f{x,y,z) = 0, as such points yield 
spatial singularity (the singularity at the t — constant hypersurface) . 
If such points exist, these can not be ignored as ignoring these will 
lead to the physical inconsistency. So, absorption of / in coordinates 
{x,y,z) is possible, only when (i)/ 7^ or (ii) / < 00 at all points 
of the t = constant hypersurface. It means that this question can be 
decided after getting the explicit definition of the function f{x,y,z), 
which is obtained in what follows. 

Using the geometrical aspect of the Ricci scalar R, R is obtained for 
the line-element, given by eq.(6.1.22), as 



dS^ = dt^ - a^{t)f{x, y, z)[dx^ + dy^ + dz% 



(6.1.22) 



f{x, y, z)[dx'^ + dy'^ + dz^] = dX"^ + dY'^ + dZ"^. 



(6.1.23) 
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(6.1.24a) 



with 



<R>=67] 



+ 



a) r 



(6.1.246) 



where < i? > is the homogeneous part of R, being its vacuum expec- 
tation value. 

When R acquires states < R >= 0, in the symmetric state at tem- 
perature T greater than or around Tc as well as < ^ >= (l/2)T'c at 
T <^Tc {as mentioned above ), the equation (6.1.24) imphes that 



^ ^ ^ ^ ^dp_Y ^dp_Y /df_Y^ ^ 

\ dx^ dy^ dz'^ I dx ) \ dy ) \ dz ) i 

(6.1.25) 

The exact solution, of the partial differential equation (6.1.25), is 
obtained as 



f = [{x - xo){y - yo){z - zo)f' . (6.1.26) 

Thus, in the state R —< R >— 0, the line-element (6.1.22) is ob- 
tained with a2(t) and p{x,y,z), given by eqs.(6.1.13) and (6.1.26) 
respectively. In the state < R >= (l/2)Tc, it is obtained with a'^{t), 
given by the equation (6.1.19) and p{x,y,z) from eq.(6.1.26). 

At this stage, it is important to test for the anisotropy and inhomo- 
geneity of the models. 

A model is isotropic, if Killing vectors exist at every point, stisfying 
equations [47] 



(6.1.27) 
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where are components of the Kilhng vector and semi-colon (;) stands 
for covariant derivatives. 

With p{x,y,z), given by eq. (6.1.26), it is found that at the point 
(xo, yo, zq) as well as in the planes x ^ xo,y ^ yo,z ^ zq, {^^.-^ + ^^-jj) 
are divergent showing singularity of these models at these points. It 
means that at these points isotropy of the model is broken leading to 
inhomogeneity also. 

Further, it is also important to see whether these singular points of 
the model are real physical singularities. This aspect can be tested 
by finding the scalar polynomials of the curvature tensor [47] . For the 

models, obtained above, it is found that 

3a 8 
a 



HI/ 



X 



49a2 
3 



[{x - Xo){y - yo){z - zo)] ^^'^ 



+ 



+ 



{x - Xq)^ (y - yof {z - ZqY 

1 1 

+ 



1 



{x - Xq) {y - yo) {y - yo) {z - Zq) {z - Zo) {x - Xo) 



3(aa + d^) + 2aa 



(6.1.28a) 

and 



R''''''R,.p. = a-'[{{x - xo) {y - yo)}-'/'{z - ^o)"''/' + •••]■ (6.1.286) 

Divergence of R^^^R^,^ and R^''^'"' R^i^pcr at above mentioned points 
shows that these points are physical singularities. 

It is interesting to see that models become homogeneous at large 
scales, when x ^ Xo,y ^ yo and z Zo- 

Second type of in homogeneous models 
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Normally, two types of models are used for manifestation of gravita- 
tional interactions in the Nature (i) static space-times representing ge- 
ometry around compact objects e.g. Swarzschild space-time, Riessner 
- Nordstrom space -time etc. and (ii) cosmological models exhibiting 
expansion of the universe. Former type models ignore expansion of 
the universe, compact objects exist and latter type of models do not 
account for the gravitational effect of these compact object treating 
these as extremely tiny particles. So, it is very natural to think of cos- 
mological models of the universe exhibiting local gravitational effect of 
compact objects as well as expansion of the universe simultaneously. 
Motivated by this idea, the second type of models are derived here [23] 
using the method demonstrated above in this section. 

Here, the line-element for cosmological models of the early universe 
is taken as [23] 



where < r < oo, < ^ < tt and < < 27r. 

In this case also, f^(r) can not be absorbed in [dr'^ + r'^{d9'^ + 9'^d(f)'^)] 
at points, where f^(r) vanishes or it is divergent. So, hke the case of 
first type of models, this question is resolved after getting explicit form 
of f^(r), which is obtained as follows. 

Here, R is obtained for the hne-element, given by eq.(6.1.29), as 



with < R>, the homogeneous part of R, given by eq.(6. 1.24b). 



dS^ = dt'' - a''{t)f\r)[dr^ + r^{d9^ + 9^d(f)% 



(6.1.29) 




(6.1.30) 
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So, in the vacuum states R =< R >= 0, at temperature T greater 
than or around Tc (in the symmetric state) as well as R =< R >~ 
(l/2)T'c at T <g Tc (when symmetry is broken ), it is obtained from the 
equation (6.1.30) that 



which is integrated to 



1 - 



To 



1 4/7 



(6.1.31) 



(6.1.32) 



Thus, in the state R —< R >— 0, the line-element (6.1.29) is ob- 
tained as 



\t\ 



1-15 



4/7 



[dr^ + r^{de'' + eH(t)'')l (6.1.33) 

when the temperature T Tc. 

In the state R —< R >— Tc/2, the line-element (6.1.29) is obtained 

as 



dS^ — dt^ — a^exp 



{t-tc 




1 - 



4/7 



[dr'' + r''{d0'' + 9^dcP% 
(6.1.34) 



when the temperature T <^ Tc. 

In these models , it is obtained that 



R' 
R^^R/j.^ 



49r6 
144r2 

49r6 



1 



1 - 



r 



-36/7 



+ 



-36/7 



+ 



(6.1.35a, 6) 

These equations exhibit divergence of scalar curvature polynomial of 
curvature tensor (which will contain R'^ and R^^'^R^„) at r = and 
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r = tq. It means that models, given by eqs.(6.1.33) and (6.1.34), possess 
real physical singularities at the point r = and the surface r = tq 
[47]. 

6.2. Creation of spinless particles in homogeneous models 

Manifestation of the Ricci scalar as a physical field also, motivates 
to study its interaction with another scalar field 0. The action for such 
a theory is 



(6.2.1) 

where ^ri~^R = with ^ being the dimensionless non-minimal cou- 
pling constant, V^{R) is given by eq.(6(i).6b) and ^XR^cjP' is the R — (f) 
interaction term with the dimensionless coupling constant A. As the 
model is homogeneous, R —< R >, which is the vacuum expectation 
value. 

The invariance of S*, under transformation — >■ -|- 50, yields the 
Klein-Gordon equation for as 

{U + in-^R + \R^)(P = {). (6.2.2) 

For the purpose of the second quantization, the general solution of 
the equation(6.2.2) can be written as 



(t> = J2[akfk{t)exp{ikaX^) + alf;{t)exp{-ikax% (6.2.3) 

k 

Particles are produced due to conformal symmetry breaking through 
the mass term . In eq.(6.2.2), the mass term {^r}~^R + XR^) vanishes 
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for the state < R >= 0. So, no production of particles are expected in 
this state. 

In the true vacuum state < R >— {l/2)Tc, the geometry of the model 
is given by the line-element(6.1.19). Connecting eqs.(6.1.19), (6.2.2) 
and (6.2.3), the differential equation for the kth mode is obtained as 



/fe = 0, 
(6.2.4) 

which can be re-written as 



with 



|e.p{-(<-g/|}+l£,-'r. + iAT|-ag^]„. = 0(6.2.5) 



fk = Vkexp{-J ^-^{t -tc)). (6.2.6) 

The in — state is obtained in the extreme past i.e., when t — > — oo 
and the out — state lies in the extreme future i.e., when t — > +00. So, 
(f) can be written in terms of in— and out — states as 

^ = J2Kfr(t>Mihxn + a^rfr{t)exp{-zk,x')^ 

k 

= J2K'''fr\t)exp{tkax')+al"'''fr'{t)exp{-tk,x% 

k 

(6.2.7a, 6) 

The normalized solution of eq. (6.2.5), in the in — state, is obtained 

as 

t;^ = exp{+iuk{t - tc))/y/2aluk (6.2.8a) 

with 



nu 



^ = ^ert + ^Ar,^-^ (6.2.86) 
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It yields 



l27Tr, 



C = exp{+i[uk - ^](^ - Q + ikaX^))/V^G^^^k (6.2.8c) 

connecting eqs. (6.2.6), (6.2.7) and (6.2.8a). 

Similarly, in the out — state, the normalized solution of eq. (6.2.5) 

is obtained as 



= exp{-iuk{t - tc))/y/24^k (6.2.9a) 
with i/jfc given by eq. (6.2.8b), which yields 



/ 27T 

= exp{-i[uk - ^J^Kt - tc) + ikaX^))/V^6^^^k (6.2.96) 

on connecting eqs. (6. 2. 6), (6.2.7) and (6.2.8a). 

The in— and out — states belong to the same Fock - space, 0^"* can 
be represented as linear combination of 0^" and 0^'" as well as 0^" can 
be written as linear combination of 0^"* and Hence [34], 

(/)r* = + /3fc0r, (6.2.10a) 
where Bogolubov coefficients and (3^ are given as 

Jt=tc 

(6.2.106) 

and 



Pk = (0r , C") = ^ / d'xa%t)[<f>r\dtcl>T) - {dtcprm (6.2).10c) 

Jt=tc 

employing the scalar product (0i, ^2) defined as 

(01,02) = -W d^X^/-gt=constant[(|>l{^t(|>*2) " (5t0l)02] 

</ t=constant 
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with 

(c,c) = i,(0r>r) = -i- 

Moreover, Bogolubov coefficients satisfy the relationship 

\ak\^-W^l. (6.2.11) 

Using eqs.(6.2.8c) and (6.2.9b) in (6.2.10b) and (6.2.10c), it is ob- 
tained that 

— and = 1. (6.2.12a, 6) 

Eqs.(6.2.11) and (6.2.12) imply 

vl = ^ (6.2.13a) 

on using eqs.(6.2.12) in eq.(6.2.11). Further, connecting eqs. (6.2.8b) 
and (6.2.13a), the result is obtained as 

r,= (12^. (6.2.13^ 

The relative probability of production of a particle - antiparticle pair, 
in mode k is given by 

employing eq.(6.2.13b). 

This result is true for all modes, being independent of k. So, the 
absolute probability of creation 1, 2, 3, • • • particle - antiparticle pairs 
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is given by [50, 51] 



oo oo 



J2 lout < 0\n\0 >in \^ ^ J2 H ^ 
n=0 n=l fc=— oo 

oo oo 



n E 2(n+l) 

fc=— oo n=0 



(6.2.15) 

on using eqs.(6.2.12a),(6.2.13a) and (6.2.14). 

It means that production of infinitely many particles, in all modes, 
is certain. The energy density of these particles can be caculated as 

f 0° =^n< 0|T°| -out < 0|T°| (6.2.16a) 

where 

= d''(t)*d^(j) - hd''(l)*di,(l) - iCv'^R + XE?)(j)*(l). (6.2.166) 

Anamolous terms get cancelled in eq. (6.2. 16a), as Tq is the difference 
of Tq in in— and out— states . 

6.3. Creation of spin-1/2 peirticles in homogeneous models 

Due to the physical nature of the Ricci scalar, its interaction with 
spin-1/2 particles is also important. 

Jd'xV^ [i(7^'^a^i?a,i?-y^(i?)+^(V'^7'^VM^-5W+c.c.)], 

(6.3.1) 

where gRipip is the R — ip interaction term with the dimensionless 
couphng constant g and V'^{R) is given by eq.(6.1.6b). As the model 
is homogeneous, R —< R >, which is the vacuum expectation value. 
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Here = — T^. Dirac matrices 7^ as well as are defind in 
eqs.(4.8), (4.9a,b), (4.10) and (4.10a,b). 

Invariance of S (given by eq.(6.3.1)) under transformation 'ip ^ 'ip + 
5'ip yields the Dirac equation 

(i^''D^-Mf^i; = Q, (6.3.2a) 

where 

Mf = gR. (6.3.26) 

(a) Solution of the Dirac equation 

In general, solution of the Dirac equation i/j can be written as 

'^^J2Y1 {hk,s)Aik,s) + , (6.3.3a) 

s=±l k 

{^(ik,s)^\,s) + ^{nk,s)^^d{k,s)) , (6.3.36) 

s=±l k 

where 

^{ik,s) = f{k,s){t)e~'^''''^Us (6.3.3c) 

and 

il^{iik,s) = g{k,s)it)e'^''''''us. (6.3.3(i) 

with g(^k,s){t) = f{-k,s){t) 
In eqs.(6(iii).8c,d) 

u( = (1000), = (0100) 

= (0010), u^^ = (0001), (6.3.4a, b, c, d) 

where the index (T) stands for transpose of the column matrices u±s 
and u±s. In eqs.(6.3.3a,b), 6it,s(^I,s) are creation (annihilation) opera- 
tors for the positive - energy particles and d-k,s{d^-k,^ are creation (an- 
nihilation) operators for the negative - energy particles (anti-particles). 
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The geometry of the model, under consideration is given by the dis- 
tance function, given by eq.(6.1.9) with a^(t) defined by the equation 
(6.1.19). The equation (6.1.9) is re-written in terms of conformal time 
as 

dS'' = (f ) [df'' - dx'' - dy^ - dz^] , (6.3.5a) 

where 

/"* dt' 

f= — — and A{f)^a{t). (6.3.56, c) 

J (^{t ) 

Connecting eqs. (6.3.4) and (6.3.5), tetrad components are obtained 

as 

e° = e} = = eH ^4^- (6-3.6) 
As a result, Dirac matrices 7^* are written as 

r = :^^rK (6-3.7a) 

and the Dirac equation(6.3.2a) reduces to 

[i{fdo + fdi + fd2 + fd^ - gA{f)R]iP = 0. (6.3.76) 

Employing eqs. (6.3.3) to the equation (6. 3. 7b), following equations 
are obtained 



[i{fdo + Tda) - ^gA{f)T,],p^ik,s) = 0. (6.3.8a) 



[tifdo + Tda) - ^gA{f)T,]^(^jik,s) = 0. (6.3.86) 
for mode k and spin s in the state < R >= (l/2)Tc, when T <^ Tp. 

Operation of [—i{^^do + ^°'da) — ^gA{f)Tc] on equations (6. 3. 8) yields 
ordinary differential equations 



^ + + 'fndoAif) + ^A^(f))/(M = (6.3).9a) 
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9{k,s) 
drf 



+ (^k^ + ^TeaoA(f) + f^A\f))g^k,s) = (6.3.96) 



using 



7°Ms = eUsi j^Us = and e = ±1. (6.3.9c, d, e) 

As g{k,s) — f{-k,s), it is better to concentrate on one of the two 
equations(6.3.9a,b). So eq. (6.3.9a) is chosen for further analysis. Con- 
necting eqs.(6.1.19) and (6.3.7) 



dt' 



-exp 



-{t-tc 



247]. 



^[-H{t-t,)] 

aM 



with = Tj2Ar] and 



A{f) = a{t) 



Ht 



(6.3.10a) 



(6.3.106) 



Using A{t), from eq.(6.3.10b), eq.(6.3.9a) is re-written as 



(Pf{k,s) 
df2 



2rr2 



+ (*^ + i^ + ^)/<M-0 (6.3.11) 
Eq.(6.3.10a)imphes that the open interval —oo < t < oo corresponds 
to the open interval — oo < r < 0. The equation(6.3.16) yields solutions 



(6.3.12a) 



= re^'''^^2ikf F^{1 -v-\- 2v- ^2ikf) 

where 



(6.3.126) 



u — - 
2 

1 
2 



, ieoTc g^T^ \ ^ 1/2 
1 ± U - 4 + — — ' 



(6.3.12c) 
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as gTc » H. Here iFi{a\h\x) is the confluent hypergeometric func- 
tion, wliicli is connected to tlie Bessel's function tlirougli an identity 
[48] 

Using tlie identity (6.3.13) in eqs.(6.3.12), tlie solutions are obtained 

as 



J{k,+i) = (7r/4) ' e"- A / y >^ i J^^^gT,/(2H+i/2){kr) 

(6.3.14a) 

= (7r/4)V2e(^-^)(V4)±...T./(4..)^i/2^^^^^^^^^^^^^^^(^-) 

(6.3). 146) 

Using the asymptotic form of the Bessel's functions [49] 



(6.3.15a) 



it is obtained that 



_ J(2fc)-^/2g-^*^^ forA;>0; 
ri-oo •^("'+^) " \ (2A;)-V2e+^*^- for A; < 0, 

which are normalized. It shows that the solution(6(iii).19a) is nor- 
malized for f — oo. The wronskian of this solution(6.3.14a) is con- 
stant. It implies that the solution(6.3.14a)is not only normalized for 
f — > — oo, but for all f in the interval (— oo,0]. The same result is 
obtained iorf(k _iy 
For small value of z, 
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Ji{^) ^ (n^- (6(m).206) 
For analyzing asymptotics of solutions for t — > ±00, eqs.(6.3.14) are 
written as 



_(9^+i)/2l 



M+i) - <| A,f'/^J_is^^^^,^{kf) as t ^ -00, 

(6.3.16a) 



are written as 
/(fc-i) 



A2T^/^J(9^_(_.)/2(^''') as i ^ 00 
>l2T^/^J_(92;+i)/2(^^) as i ^ -00, 

(6.3.166) 

Also using g^k,s) = f(-k,s): 

{^if^/^J(9ii_^.)/2(^^^) as i ^ 00 
AifV2j_[k^.)/2(-^^) as t ^ -00, 

(6.3.16c) 

A2r^/^J(9^^.y2(~^^) as i ^ 00 
>i2T^/^J_(9^^j^/2(~^^) as i ^ -00, 

(6.3.16d) 



9(k,-i) 



In eqs.(6.3.16a,b,c,d), 

^ Qy/'e^^^-l'^^^l/^) (6.3.16e) 

and 

^2= (6.3.16/) 
Thus the normalized solutions '4'{ik,s) and ip{iik,s) are obtained as 



^(^M) = ^ (^) (6.3.17a) 

and 
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Arrk,s) = ^ (^) '^%ik,s)e(^''^'^us (6.3.176) 
in the three-space volume V, using the normahzing prescription [34] 

{'4^{k,s),'^{k',s')) = / d^X\/ -gt=constant'4^{k,s)l^'^{k' ,s') = ^ss'^kk'- 

J t=constant 

(6.3.17c) 

(b) Production of particles 

For the study of creation of spin-1/2 particles, it is better to go 
to Fock-space formulation, where the m-state {t — > — oo) as well as 
out-siaXe {t — > — oo) is defined as 



blls)\0 >in= V[-k,s)\0 >in^ and i„ < 0|0 >i„= 1. (6.3.18a, 6, c) 
and 



(6.3.19a, 6, c) 

The decomposed form of ip is written in the in-region as well as 

OMt-rcgion as 



^ = EE(^(MV'(/M) + (^r-^,-.))Vr.(-.,-.) 

s=±l k 

E E(^(M^°"m) + ((^(-fe,-s))V//(-fe,-s)) 



s=±l k 



(6.3.20a, 6) 



and 
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s=±l k 
s=±l k 

(6.3.20c, d) 

Bogolubov transformations for Fermi-fields are written as [52] 

d{-k,-s) = <^(A;,s)(^(fc,s))^ + /^(A;,s)(^(fc,s))^ 

(6.3.21a, 6, c,d) 

Bogolubov coefficients a(^k,s) and satisfy the condition [48] 

|«(mI'+IAm)I' = 1, (6.3.22a) 

which are given as 

«(M= / d'x^-gt=constanti^t^k,sMrls)y (6-3.226) 

•/ t=constant 

and 

A^,.) = / C^'^V-^t=con.tantV'}7(-fc,-.)(V^77(-fe,-.))^- (6.3.22c) 

•/ t=constant 

Connecting eqs.(6.3.21), (6.3.22) and (6.3.27), 

(7r/4)eW2[_((24ry)rjV-^e-^ft^-^-)]^- fky^ 

«(M - [1 _ i(i + 5re///)]v^[i + i(i + 9TJH)Y/^ ^'-'-''"^ 



^(M - [1 _ 1 (i + ^T.///)]V2[i + i(i + gTjH)]y- I 2 J ^'-'-'''^ 
Eqs. (6.3.23) yield 
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, ,2 _ (7r/4)^e-[{(24r^)TjV^e-^fa-*-)]^ 

' ([(1 - gVWc)' + i/4][(i + + i/4]V4 ^^■^■'^"^ 



([(1 - ^V6^)^ + 1/4][(1 + ^V6P;)^ + 1/4] V4 ^'•'•'^''^ 
The absolute probability of no creation of particles are given as [50, 
51] 

-2 



\out < 0|0 >in r = n 

k 

= ea;p[^/n(|a(fc,+i)|-V(fe-i)r^) 

(7r/4)4[{(24r/)Te}i/4e-^(*i-*=)]^ 



([(1 - gvWcr + 1/4] [(1 + gvWcr + 1/4]^^ 

(7r/4)^[(24?7)7;e-^^(*i-*°)] 
~ (7r/4)^(4//)e-^^(*i-*=). 

(6.3.25a) 

The decay of the in-state per unit time per unit volume is given by 

f = -ViHn\out < 0|0 >i„ I' 

~ a;3e[-3^(*-*-)][4i/(t - g - /n(7r/4)^(4//)]. 

(6.3.256) 

It shows that when t > tc, the decay of the in-state is very fast. As a 
result, particle- antiparticle pairs will get created. 



6.4. Creation of spinless peirticles in inhomogeneous models 
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Here, inhomogeneous models, given by distance functions (6. .33) and 
(6.1.34), are considered for states < R >= and < R >= Tc/2 respec- 
tively. 

Action for the scalar field is taken as [53] 



1 



The Klein - Gordon equation, obtained from this action is 



(6.4.1a) 



[□ + M2]$ = 0, 



(6.4.16) 



where 



= A<R>^ 



(6.4.2) 

In the background geometry of the hue - element (6.1.33), the equa- 
tion (6.4.1b) is re- written as 



3d 9$ 1 r ro 
dt'^ a dt a^l r J 



-4/7^2^ 1 



Qj.2 j.^^^ 



--] 
r J 



1 d / . a 



To 



-6/7 d 

dr 



2/7- 



sinO do 



+ 



r J i dr 



$ + M^^ = 0, 



(6.4.3) 

where dot (.) over the variable denotes the derivative with respect to 
time. 

The solution for this equation is taken as 

3/7-1 -1, 

k,l,m 



"(27r)i/V(l-^)' ]" J][Afe,^*,,„(t)ll„(^^,^/.)e^'=-'- + c.c." 

k,l,m 

(6.4.4) 

where r > ro,m — —I •■ - + 1,1 — 1,2,3, •■ ■ and — 00 < k < 00. Yim{9, ip) 
satisfies the equation 
1 d 



sinO 89 



d\ Id' 

+ 



Yim^-l{l + l)Yim (6.4.5a) 
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with normalization as 



(6.4.56) 



Connecting eq. (6.4.3) with eq. (6.4.4) and using eq. (6.4.5), it is ob- 
tained that 



Ak.r 



T "^^ T 

a 



A-2 ^ /(/ + 1) ^ 6ro(2r - tq) 



q2 q2j,2 Ya'^r^(^r — vq) 



= 0. 



1 



To 



-4/7 



(6.4.6) 



Using the convolution theorem [54] and integrating over r, eq. (6.4.6) 
looks like 



f ikr, \-i 3d • r r° 



V Jro+e 



1-^ 



-4/7 



+ 



/■oo ^-ik.y 




Jro+e 


L y J 



-4/7 , , 6ro 
dy+ 



- /■(-)) L r 



7a'^vJro+ey^iy -^o)^^ y 



= 0. 



-4/7 

dy 
(6.4.7) 



Here e is an extremely small positive real number. 

Details of the evaluation of integrals with respect to y are given in 
the Appendix B. Using these results, in eq. (6.4.7), one obtains 



T '^'^ T 

^ klm ^ Mm + 

a 



= 0, 



(6.4.8a) 



where 



Xkim = {l/a'rj)e'/\ro + e^/'cosikiro + e)} 



(ro + ey 

(6.4.86) 



Case 1: The case of state < R >— 



Using a{t) from the line-element (6.1.33) for this state, the equation 
(6.4.8a) is written as 
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which yields the solution for t > as 



-1/2 



(6.4.9) 



(6.4.10a) 



AJ_i/2{t) + BY_i/2{t) 

where A and B are integration constants. Here r is defined as 

T = ^JtpXk^m(t + tpajo) (6.4.106) 

showing that r — > cxd when t — >• oo. 
For large r, 

cos(r) sinfr) 
J_i/2(t) ~ and 1^i/2(t) ~ 



VW2 ' a/W2' 

So, when r is large, 

= [7TT^/2]-^^^[AcosT + Bsinr] 

[1 + i)e-'^ + {1 - i)e'^ 



\ 1 




-2tpXkijn. 





using the normalization condition 

(*«m,*«m) = 1 = -($«„,$« J, 

where 

ro\ 1/71-1 



(6.4.11) 



(6.4.12a) 



(27r77-ipF)^/V(l-^) ]" *«„(t)e^'=-'-ll^ (6.4.126) 



V^dEi^klmidt^U'm') 



and the scalar product is defined as 

{^klm, ^klm) ^-i II 
Jro+eJo Jo 

-{dt^klmWk'l'm'Wlmie,^', {6A.12c) 

where E is the t — ti hypersurface, ^/—ffE — 
sinOdrdOdip. 



6/7 

1-^1 and = 
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Thus 



^out _ 
klm 



1/2 



.2tpXkim 

For t < 0, eq. (6.4.9) yields the solution 



T-\l + i)e- 



(6.4.13) 



(6.4.14a) 



A'j_i/2(f)+5'F_i/2(f) 

where A' and B' are integration constants. Here f is defined as 

(6.4.146) 



— —\JtpXum{—t + tpal) 



showing that f — > — oo when i — > — oo. 

Using above approximations for Bessel's functions and normalization 
prescription, it is obtained that for t — —oo 



klm 



klm 



1/2 



+if 



(6.4.15) 



Since ^^/i^ and ^^^im both belong to the same Hilbert space, so one 
write 

- E [^°.?m*»m(^, V') + ^.m(^, V')^ . 

k,l,m 

As a result, one obtains 

= + Pkim^tL, (6.4.16a) 

where 0;^/^ and (5kim are Bogoliubov coefficients satisfying the condition 
[31, 34, 52] 

|2 



\'^klr. 



1. 



The in- and out- vacuum states are defined as 



4in 

-^klm 



Moreover, 



in >=0 = A-^|out>. 



/lout A\n , o A- 

^klm — "'klm^klm Pklm^. 



l*in 

^klmi 



(6.4.166) 



(6.4.17a, 6) 



(6.4.17c) 
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Connecting eqs.(6.4.11) - (6.4.17c), it is obtained that 



Pklm — 2^ 



-i(r+f) 



-i(r+f) 



\— + alQtpJ \ ti ■ 



no 

,2 



These results yield 



implying 



5ti 



(6.4.18a, b) 
(6.4.19a) 
(6.4.196) 



Eq.(6.4.19b) shows that when ti is sufficiently larger than 5tp, \/3kim\'^ — 
0. It means that spinless particles will be created in the state < R >— 
only when ti < 5tp. For ti sufficiently larger than 5tp, there will be 
no production of scalar particles, in this state. 

Anamolous terms emerge when vacuum expectation value of the 
energy-momentum tensor components is evaluated for if dimen- 
sion of the space-time is an even integer. These terms get cancelled for 
the difference of vacuum expectation values T^j, in in- and out-states. 
Thus T^jy for created particles is defined as [51, 53] 



< 7),^ >=< out|r„,,Jout > - < in|r„,,Jin > . 



ixu I 



(6.4.20a) 



T^i, for (f) are obtained from the action ( 6.4.1a) as 



(6.4.206) 
with its trace 

T = -(1 - 12riAR)d''^*d^^ - UrjA.^ < R >^ (6.4.20c) 
Now trace of T^j,^ for created partcle-antiparticle pairs is obtained as 
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< T >=< out|r|out > - < in|T|in > (6.4.20d) 

Eqs.( 6.4.12b), ( 6.4.13), ( 6.4.15) and ( 6.4.20 d) lead to trace of 
created particles in the state < ^ >= as 



< 7\(6) 0. 



(6.4.21) 



Case 2: Tlie c-ase of state < R >= 



Using a{t) from the line-element ( 6.1.34) for this state, the equation 
( 6.4.8) is written as 



T 3 /Te- 



klm ^-(t-to)^/Tc/6v 



^Mm = 0, (6.4.22) 



where — \A.T^ using the definition of Mh from eq.( 6(iv).2). The 
equation ( 6.4.22) yields the solution for i > as 



klm 



TT 


1/2 


"67/" 


sinh2Tra. 




T 



J±2ia (7A:^me-^(*-*-) V^) , (6.4.23a) 



where 



and 



P 677 rl,. 3Tci 
-4 ^ 32r] 



2 _ 2AriXkim 

For t < 0, the equation ( 6.4.22) yields the solution as 

na 



(6.4.236) 



(6.4.23c) 



klm 



- 1/2 


67]" 




T 

L J. g J 



1/4 3, 
g4^ 



(t-tio)V^c/6?7^ 

JT2.a(7«rne^(*-*-)V^). (6.4.24) 
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Using the approximation of 

"^"^^^ - 2"r(l + n) 
for small x in eqns.( 6.4.23) and ( 6.4.24) , for t ^ oo 
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^out 
klm 



rjsinh{2Tia' 



■ 1/2 






T 



1/4 



r(l±2m)/V 2 



±2ia 



X 



-{t-tw){^T^a)^/Tc/QrJ 



(6.4.25a) 



and for t — oo 



— 

klm 





1/2 


'6?]' 


r)sinh2Tra. 




T 



\r(1^2ia))\ir) ^ 

(6.4.256) 



= (t-tio)(|=Fi«)\/^c/6r, 



Using ^^J'^ and one obtains 



(Xklr 



and 



implying that 



3i 
9 



'feim 



(6.4.26a) 
(6.4.266) 
(6.4.26c) 



which ensures creation of spinless particles. Connecting eqs.(6.4.16c) 
and (6.4.26b) 



12 13 

\Oiklm\ — 



(6.4.26d) 

Moreover, the absolute probability for no particle creation is given as 



< out in > 



1-2 



klm 



9 + 16a2 ■ 



(6.4.27) 



Eqs. (6.4.26 d) and (6.4.27) yield 



= 1/4. 



(6.4.28) 
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Connecting eqs. (6.4.2), (6.4.23b) and (6.4.28), it is obtained that 



Ml 



16 \ Qr] 

Using the convolution theorem, eq. (6.4. 12b) can be written as 



(6.4.29) 



{2pi) 



-1/2 -ikr 



{2'pi) 



-1/2 -ikr _ 



.4/7 



y{i - '^)Y\^''ydy\-^kUt)YUeA) 



7 

cos[k{rQ + e)]^ kim{t)Yim{() 1 0)- 



(ro + e)iVv 

(6.4.30) 

Connecting eqs. (6.4.20 d), (6.4.25) and (6.4.30) as well as taking 
average over 6 and 0, trace of the energy-momentum tensor for created 
particles, in the state R — (l/2)Tc, is obtained t — t2e as 



r2(,) = J](l/27r)(|)'^m/.[(3/4)^/^(t2e-W] 



klm 
X 



(ro + e)2Vv^^^^(^o + ^) [(13/16)^(6r;AT, - 1) 



(6.4.31) 



6.5. Creation of spin-1/2 pcirticles in inhomogeneous models 

In the homogeneous model (6.1.33) and (6.1.34), the general solution 
of the Dirac equation (6.3.7a) is taken as and -0^, given by eqs. (6.3.8a) 
and (6.3.8b) with 



Ik.s 



(27r77)i/V(l - ^)'^']"'5^A;^,.(t)ll^(^,0)e'^-X (6.5.1a) 

l.m 
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and 

l,m 

(6.5.16) 

with gklm,s{t) = f-klm,s{'t)- 

In eqs.(6.5.1a, b), Ug and its are given by eqs. (6.3.9a, b,c,d). Moreover, 
m — —I, • • • , +1; I — 1,2,3, ■ ■ ■ and — oo < k < oo. 

Connecting eqs. (6.3.7a) and (6.5. la,b), it is obtained that 

(iYD^-Mf'^tl;ik,s^O, {6{v).2a) 
Now using the operator ^ — i^^D^ — Mf^ from the left in eq. (6.5.2a) 



as 



it is obtained that [23] 

(n + (4r?)-^ < ^ > < ^ >2 ^i)ik,s = 0, (6.5.3a) 

Similarly, from eq. (6.5.2b), it is obtained that 

(u + (4r^)-i <R> <R>^ )iJiik,s = 0. (6.5.36) 

When temperature T is not very much below Tc one obtains the 
vacuum state < R >— 0, where eqs.(6.5.9a,b) reduce to 

0^ik,s = = □V'77fe,. (6.5.4a, 6) 

Particles are created due to conformal symmetry breaking, which is 
caused by the mass term. In the state < R >— 0, there is no term 
to break this symmetry. So, production of spin- 1/2 particles are not 
possible in this state. 
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As given in section 6.1, one obtains the state < R>= ^Tc when the 
temperature falls sufficiently below Tc. In this state, the universe obeys 
the geometry given by the line-element (6.1.34) . 

In what follows, creation of spin- 1/2 particles is investigated in the 
model given by this line-element. In the background geometry of this 
model, for every k, 1,171 and s, eq.(6.5.3a,b) look like 



„ik.r 



jklm.s ' J klm,s ' I 1 

a Via 



^ /(/ + !) ^ 6ro(2r-ro) 
^ a^r^ 7a^r^(r — ro)^ ■ 



1- 



-4/7 



= 



r 

(6.5.5a) 



connecting eqs.(6.5.1a) and (6.5.3a).. Here 

M| = (477)-^ <R>+g<R>' 



o . ■ (6.5.56) 

8r] A ^ ^ 

Using the convolution theorem [54] as above and integrating over r, 
one obtains that 



poo 
Jro+e 



3d 



fklm,s ~l~ fklm,s ~l~ ^ o 



-ik.y 



^ VJro+e 



ro+e 



y 



-4/7 



6rn 



-4/7 



dy 



(2y - rp) 

y'^iy - rof 



y^ 



-4/7 



0. 



(6.5.6) 



e-"^-ydy + M]jfki„,. 
Details of the evaluation of integrals with respect to y are given in 
the Appendix B. Using these results, in eq. (6.5.7), one obtains 



J klra.s "I" 1 / J klra.s ~r 

2 y 677 



«?0 



where 



Xkim,s = {l/a^n)€'l\To + eY"Cos{k{To + e)} 



fkim,s = 0, (6.5.7a) 
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-yroe-'(ro + e)-^ + ^l^-^r^ + e)" 



(6.5.76) 



Eqs.(6.5.8a,b) yield the solution for i > 0, 

him,s - Cie±t(*-*-)x/^ J±2.« (7...,.e-^(*-*-)x/^) , (6.5.8a) 
where Ci is a normalization constant, 



and 



a 
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M' 



2 3Te 



6^rl 2^.2 , 5T, 



^ 3277 J 

using the definition of M| from eq. (6.5.5b). 
Connecting eqs.(6.3.9a,b,c,d) and (6.5.8a,b) 



(6.5.86) 



(6(v).8c) 



Ik.s 



Lm 



X 



J±2ia (7«m,.e-^(*-*-) 1^^(^, 4>)e'^-us. (6.5.9) 

Normalization of ■07fe,s is done using the prescription (6.3.22c) at the 
t = tc hypersurface denoted as E onwards. As a result, Ci is obtained 

as 



J±2ia 



(6.5.10) 



Thus 



Z,m 



klm,s J 



Similarly 

l,m 



g-!(*-tc)V^V6^^X 

(6.5.11) 



^IIk,s = [(27r77)^/V(l-^)'''y5^|j±2^a(7- 



klm.s I 



-1 



g-|(*-*c)v^V6^x 
(6.5.12) 
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Asymptotics of these solutions, when t ^ oo, yield 



X 



l,m 



r(l - 2ia) 



e+i'-it-tc)V^Yim{e, (t>)e'^-'us. (6.5.13a) 



and 

,/,out 



(6.5.136) 



For t — > — oo, one obtains 



(27r77)^/V 1 - ^ Y,\j±2ia(lklm,s) 



3 
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X 



r(l + 2^a)l 72:^^ ''"e-^"(*-*=)v^y;^(^»e^'=-X. (6.5.14a) 



and 



(27r77)VV(l - ^)'Y'Ek±2^«(^-^'-.^)r'^^^*"*^^^^x 



l,m 



r(i + 2m) 



-2ia 



^(t-t.)V^YU9,<f>)e-''-'-u,. 

(6.5.146) 

In eqs.(6.5.13a,b) - (6.5.14a,b), r(x) stands for the gamma function of 

X. 

The decomposed form of t/j and ^'^ are given by eqs.(6.3.25a,b) and 
eqs.(6.3.25c,d) respectively. 

Connecting eqs.(6.5.13a), (6.5.14a) and (6.3.27b) 



ak.s = (1/2) [r(l - 2ia) 
sinh{27ra 



■J±2ia{lklm,s) 



r(l - 2ia) 



(6.5.15a) 
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I3k.s = (l/2)[r(l + 2ia) 
sinh{2TiQ.) 



1 -2 



V 2 



-Aia 



J±2ia{lklm,s) 



r(l + 2ia) 



(6.5.156) 

Conditions (6.3.27a) and eqs.(6.5.15a,b) yield an useful result 

(6.5.16) 
-k 



1 

4 



suggesting that a should be very small, using the symmetry k - 
caused by symmetry of the model under r —r. 

The relative probability of creation of a particle - antiparticle pair is 
given as 



^klm.s 



= 1. 



(6.5.17) 



Absolute probability of the creation of particle - antiparticle pairs 
requires the total probability of creating 0, 1, 2, • • • pairs to be unity 
[7], which means that 

Nklm,s{^ + ^klm,s + ^klm,s ) 

where Nkim,s is the probability of creation of no pair of particle-antiparticle. 
Using eq.(6(v).17), it is obtained that 

1 



klm,s 



0. 



1 + 1 + 1 + -- 

It shows that probability of vacuum to remain vacuum is 



< out in > 



klm,s 



0. 



(6.5.18) 



klm.s 



implying certainity of creation of particle-antiparticle pairs. 
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The action S, given by eq.(6.3.1), yields components of energy - 
momentum tensor for ip field as [52, 55] 
1 6S 

-ar]{Df,D^ - g^^a)'4;ijj + c.c. 

(6.5.19) 

Connecting eqs. (6.5.20a) and (6.5.19), energy density of created par- 
ticles in the state R = l/2Tc is obtained as 

p =< > = < in| - iV^7°9o - a-qRli)^) + 2aRi)'il) + c.c. |in > 

— < out| — iijj'y^do — ar]RQip'^ + 2aR^pip + c.c. |out > 
= < in| -iVS7°9o + l/2(7TcVSV' + c.c. |in > 

- < out| - i'07°9o + l/2crTc-(/'-(/' + c.c. |out > 

(6.5.20) 

using i?Q = 3a/a = Tc/2rj. 

In- and out- t/^, given by eqs. (6.5. 13a,b) and (6.5.141,b), recast eq.(6.5.20) 

as 

So; e^/"^ 

P = ll — . ^^N22/7 ^^^'[^(^o + e)]sinh[{2>/2){t - Q^/tJH 

klm ^ ^' 

(6.5.21) 

using convolution theorem for il) as given in the preceding section. 
Similarly trace of stress - tensor for created spin-1/2 particles is given 

as 



(6.5.22) 
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using In- and out- solutions iorip, given by cqs.(6(v).13a,b) and (6.5.14a,b). 
Eq.(6(v).16) shows number density of created spin-1/2 particles as 



1 



(6.5.23) 



2 



yielding energy density of created particles also as 



P = \Mf. 



(6.5.24) 



Here summation is done with the help of the Riemann zeta function 
defined as ^(s) = ^^i^t-""*- Moreover, 



klm klm 

klm kl 

= E[2C(-l) + C(0)] = -(4/3)C(0) = 2/3 



as ({—2m) = 0. Also ({0) = — | and C(~l) = ~J2 obtained through 
the analytic continuation. 

7. Contribution of riccion to dark energy and dark matter 

of the universe 

Experimental probes, like luminosity measurements of high redshift 
supernova [56], anisotropies in cosmic microwave background [57] and 
observational gravitational clustering [58], strongly indicate late time 
accelerated expansion of the universe. Theoretically, accelerated expan- 
sion of the universe can be obtained cither by modifying left hand side 
of Einstein's equations, like introduction of C-field in steady— state 
theory adhered to the Perfect Cosmological Principle [59] or by us- 
ing energy momentum tensor,having dominance of exotic matter, with 



k 
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negative pressure violating the strong energy condition. This kind of 
matter is known as dark energy (vacuum energy), which has drawn 
much attention of cosmologists today. Past few years have witnessed 
concerted efforts to propose different dark energy models. Important 
efforts, in this direction, are scalar field models like (i) quintessence [60], 
(ii) k-essencc [61], (iii) tachyon scalar fields [62, 63, 64, 65, 66, 67] and 
models based on quantum particle production, Chaplygin gas [68, 69]. 
In these models, a scalar field acts as a source of are dark energy and 
plays crucial role in the dynamical universe. But these models have 
no answer to the question "Where are these scalars coming from?" It 
is like Higg's fields in GUTs as well as infiaton in infiationary models 
of the early universe. In the cosmology probed here, no scalar field is 
required to be incorporated, from outside the theory, to discuss cosmic 
dark energy. Following arguments in [69, 70], here also, vacuum energy 
is recognized as DE. 

All these dark energy models (mentioned above) try to explain rolling 
down of DE density from a very high value in the early universe to an 
extremely small value in the current universe as it is suggested by as- 
tronomical observations. But its initial value, in the early universe, 
is different in different contexts. For example, it is ~ lO^^GeV^ at 
Planck scale, ~ at GUT phase transition and quantum chro- 

modynamics yields its value ~ 10~^GeV^. In what follows, one-loop 
renormalization of riccion [72, 73] contributes initial value equal to 
lO^GeV^ for slowly varying time-dependent DE density to the dynami- 
cal universe beginning with a phase transition at the electro-weak scale 
Mew = lOOGeV. 
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Here, there is no role of any other scalar field in the cosmologi- 
cal dynamics except riccion which emerges from higher- dimensional 
high-derivative gravity. Thus, unlike other models, dark energy is con- 
tributed by the gravitational sector. Geometry of the space-time is 
considered to be (4 -|- £)) -dimensional with topology (8) S^, where 
is a D-dimensional sphere being the hidden extra-dimensional com- 
pact space. Further, it is obtained that D ^ 6 which makes effec- 
tive dimension of the space-time equal to 10. The observable uni- 
verse is a 4-dimensional hypersurface of the higher-dimensional world. 
So far sharpest experiments could probe gravity upto 0.1 mm. and, 
in the matter sector, probe could be possible upto clectroweak scale 
Mew =^ lOOGeV, having wavelength as short as ~ 1.97 x 10~^^cm. 
It means that it is difficult to realize 4-dimensional gravity for length 
scales less than 0.1 mm. So, it is reasonable to think higher-dimensional 
gravity for length scales smaller than this scale. This idea is parallel 
to the brane-world gravity, where it is assumed that gravity is stronger 
in higher- dimensional space-time, called bulk and only a small part of 
it is realized in the observable universe [1]. 

Thus, in the proposed cosmological scenario, our dynamical uni- 
verse begins with the initial dark energy density (vacuum energy den- 
sity) pAow = lO^GeV^ and background radiation with temperature 
Tew = 78.5GeV = 9.1 x lO^'^K, contributed by riccion at Mew The 
background radiation is caused by a phase transition at Mew This 
event is recognized as big — bang. As usual, dynamical universe grows 
in the post big — bang era. In the present theory, riccion has effective 
role in the pre— big — bang era only and it remains passive in the post 
big — bang times. 
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T.l.Riccion from iV/^ ® geometry of the space-time 
Here , geometry of the space-time is given by The distance function 



dS'^ = g^^dx^dx" - l^dQ'^ (7.1.1a) 

with 

dQ^ = dej + sin^Oidel + ■■■ + sin% ■ ■ ■ sin'^e^D-i)del, (7.1. 1&) 

where gfj,i^{iJ>, = 0, 1, 2, 3) are components of the metric tensor in M'^, 
I is radius of the sphere 5"^ which is independent of coordinates x^ and 

< ^1, ^2, • • • , d{D-i) < TT and < < Stt. 
Theory begins with the gravitational action 

/ j^(2+D)^ 
d^xd^y^-g(^^+D) hQ;(4+D)-R(4+D)+7(4+D) (-R(4+D) 

— {D - 2) (^+^) (-D+4)^J' (/.i.^aj 
where G'(44.£i) = M~^^+^'' (M being the mass scale ),a(iJ^D) = oiV^^ , 
7(4+D) = 3!(D-2) Rd = Here Vd, being the volume of 

5*^ ,is given as 

= ^ T^P. (7.1.26) 

r(L) + i)/2 ^ ' 

g{A+D) is determinant of the metric tensor gMN{M, A^ = 0,l,2,---,(3 + 
D) and R(4,+d) = R + Rd- «isa dimensionless coupling constant, R is 
the Ricci scalar in and Gn — G(^4+d)/Vd- 

Invariance of S under transformations gMN — ^ gMN + ^guN yields 
[21, 22] 



— -j^Q^ {Rmn — -gMNR{4+D)) + Oi{^i+D)Hj^^ l{i+D)Hj^^ — 0, 

(7.1.3a) 
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where 



(1) ~ A 

^MN — 2-R;MAr — '^gMN^{i+D)R{A+D) — -^gMNR(i+D) + '^R{4+D)RmN, 

(7.1.36) 



and 



^MN — ^R-MN - ^9MN^{4+D)R'l4+D) " ^jj _ c^{^ {-1^9MN^(A.^-D)R\a,+D) 
+2n(4+D)-R(D+4)-RMA?" + -RfMAr} ~ -^9MN^(a,^b')^'^R\/L+D)RmN (7.1. 3c) 

Trace of these field equations is obtained as 



6(D + 3), 1 



{D + 2)M(2+ 



327r 



+ a(4+D) [2(-D + 3) □(4+D)-R(4+£)) + 2 -^^(4+L>)] 



+ T^7(4+D)(i?-2)i?f4+o)=0. 



(7.1.4) 



In the space-time described by the distance function defined in eq.(7.1.1a 



1^(4+1)) -R(4+D) = ^R — 



1 d 



-9 9 



d 



(7.1.5) 



-g dxi^ V dx'^ > 

using the definition of i?(4+D) given in eq.(7.1.2a). 

Connecting eqs.(7.1.4)- (7.1.5) as well as using i?(4+D), Q!(4+d) and 
7(4+D) from eq.(7.1.2a), one obtains in 

■(D + 2)M(2+^)yB- 



327r 



(i? + i?D) + Q;[2(D + 3)ni? 



^^D{R + Rof\ + j::{R + Rof = 0, 



12 



(7.1.6) 



which is re- written as 



[□ + ]-^R + + ^rfR^]R + n'^'d = 0, 



(7.1.7) 



where 
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m 



IGtt 



DRn 1 Q , „r, 



A 

1? 



1 



= -77 



4(D + 3)q;' 

(L> + 2)AMMyo 



IGtt 



4(D + 3) 6 



(7.1.8a, 6, c, (i) 

where o; > to avoid the ghost problem. 

A scalar field, representing a spinless particle, has unit mass di- 
mension in existing theories. R, being combination of second order 
derivative as well as squares of first order derivative of metric tensor 
components with respect to space-time coordinates, has mass dimen- 
sion 2. So, to have mass dimension like other scalar fields, eq.(7.1.7) is 
multiplied by rj and rjR is recognized as R. As a result, this equation 
looks like 

p + ^^R + m^ + ^R^]R + ^ = 0. (7.1.9) 

Using the method given above, the corresponding riccion action is 
obtained as 



-d^Rd.R 
.2 ^ 



(^j^iRR? + ]^m^R? + ^,R^ + m 



4!' 



)]}■ 

(7.1.10) 



7.2. One-loop quantum correction , renormalization of riccion 
and solution of renormalization group equations 
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Here also, one-loop quantum correction and renormlization of riccion 
is done using the operator regularization method [36] with the action 
(7.1.10). In contrast to riccion action (5.4), the action (7.1.10) contains 
another term '&R. Due to this term, in this case, aditional term and 
equation appear in renormalized lagrangian (5.8a, b) and renormaliza- 
tion conditions(5.9). Solution of renormalization group equations for 
relevant coupling constants for further use in this section are given [73] 
as 



A 

A 
1? 



m 



■''-ew ~r 



2A. 



/ 3Aewry/3_ 1 



3Aewri -1 



•dewiconstant) 

3AnwT 



1 - 



87r2 



-1/3 



1 /I _1\ 



3AewrT-i 



87r2 



(7.2.1a, b, c, d, e) 



where r = \ln{Ml^ / ji'^) with cut-off mass scale Mew Here, Aj^^ = 
Ai(T = 0) and r = at ji — Mew according to defnition of r 
given above. 

These results show that as/x^oo(T— >— oo). A— >0 and w? — > 
and \i ^ |. 

Using these limits in eq.(7.1.8a), it is obtained that 



D = Q. 

Also,eqs.(7.1.2a) and (7.2.2) imply 



(7.2.2) 
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Re = (7.2.3) 

So, from eq.(7.2.1a) 

A = Aew (7.2.4) 

at n = Mew The equation (7.2.1c) shows that i) is independent of mass 
scale /I. So, ri'^-d, being true for all fi, is obtained from eq.(7.1.8d) as 



2 
ew 



at /J, — Mew. Here 

, _ AewM^^T/eew , 10(l+450Aew) 

with 

T^iew = ^^^^ (7.2.6) 
15 

which is derived connecting eqs.(7.1.8b,d), taking the arbitrary param- 
eter Tj = /ew and putting D = 6. 

Thus r]~^'& , given by eq.(7.2.4),is an imprint of extra six-dimensional 
compact manifold ^S"^ in the 4-dimensional universe. This term has the 
dimension of energy density. Moreover, it is not generated through 
matter, but the geometry of extra-dimensional space. So, it is recog- 
nized as dark energy density p\^^ at fi = M^^, given as 



30 



'ew 



5(1 + 60Aew) 



PAew = = W- ^ew - ^ n , (7.2.7) 



7.3. Equation of state for deirk energy and Phase transition 
at the eletroweak scale 

Dark energy density pAew > obtained through renormalization of ric- 
cion at energy mass scale /i — M^-^ , is given by eq.(7.2.7). As riccion 
is emerging from geometry of the space-time without matter, can 
also be obtained as zero-point energy of riccion, given as 
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PAit = 0) = pA,„ = (27r)-3 / y/W+^47rk'dk, (7.3.1) 
where mew is the mass of riccion at the electroweak scale M^^. Eq.(7.3.1) 
shows that pAcw diverges as /c ^ oo. But pAow i^^ given by eq.(7.2.7)), 
being finite, imphes that the integral in eq.(7.3.1) should be regularized 
upto a certain cut-off mode k = k^- As so far experiments could be per- 
formed upto Mew only, so cut-off scale is taken as kc — M^-^ as above. 
Moreover, at this scale riccions heavier than M^^ can not survive, so 
"^ew ^ -^ew Now eq.(7.3.1) yields 

taking mew = Mew- The second quantization and uncertainty relation 
imply that vacuum has energy density as well as pressure [68]. Exper- 
imental probes, like la supernova and WMAP [56, 57, 58, 81] suggest 
accelerated expansion of the universe, which requires negative pressure 
for DE. So, to have consistency with recenprobes, the isotropic vacuum 
pressure is calculated as 

It yields 

(7.3.4a) 

Though in certain models, time-dependence oi cu = p/p is also pro- 
posed , but normally it is taken as a constant. So, here also, this ratio 
is considered independent of time. As a result 



uja = Pa/ Pa = PAew/PAew = -0-4- 



(7.3.46) 
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Using = Mc„ = mew and connecting cqs. (7.2. 5), (7.2. 6), (7.2. 7) 
and (7.3.2), a quardatic equation for Aew is obtained as 



(450 - ^) + [20(450 - ^) + 30 X 2818.8] Aew + 11375.3 = 0, 

(7.3.5) 

which is solved to 

Aew = -0.013352958. (7.3.6) 
Also using Mew = mew in eq. (7.3.2), it is obtained that 

PA,„ ~ lO^GeVl (7.3.7) 

In field theories, Planck scale is supposed to be a fundamental scale. 

So, it is proposed that energy mass scale fi falls from the Planck mass 
Mp = lO^^GeV to the cut-off scale Mew = lOOGeV. When it happens 
so, phase transition takes place at /i — Mew a-nd energy with density 

Pew(r) = A - Aew = 2.5 X lO^GcV^ (7.3.8a) 

is released, which is obtained connecting eqs.(7.2.1a) and (7.3.6). 

In the proposed model (PM), this event is recognized as big — bang, 
being the beginning of universe like standard model(SMU) with the 
release of background radiation having energy density Pew(r), given 
by eq. (7.3.8a). In contrast to SMU, here, released energy density at 
the epoch of big-bang is finite (it is infinite in SMU). Temperature of 
photons Tew with energy density Pew; is obtained from 



Pew(r) = ^Tet = 2.5 X lO^GcV^ (7.3.86) 



as [73] 



Tew = 78.5GeV = 9.1 x lO^^K. 



(7.3.9) 
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7. 4. Proposed cosmological scenario, Dark energy and Dark matter 



In the standard model of the big — bang theory, it is supposed that, 
around 13.7Gyrs ago, there used to be a fireball ( an extremely hot 
object), which was termed as primeval atom by Lemaitre. Our uni- 
verse came into existence, when this primeval atom burst out. This 
event is called big — bang. 

According to the proposed cosmological picture, the primeval atom 
contained riccions, being contribution of 10- dimensional higher-derivative 
gravity to the 4-dimensional world. One-loop renormalization of riccions 
and solutions of resulting group equations yield that when energy mass 
scale comes down to electroweak scale = lOOGeV, riccion con- 
tributes dark energy density pAew = lO^GeV^. Moreover, phase transi- 
tion takes place at Mg^, releasing the background radiation. This radi- 
ation thermalizes the universe upto the temperature T^^ — 9.1 x 
Here, the event of phase transition is recognized as big — bang, which 
heralds our dynamical universe having the initial temperature T^^ — 
9.1 X lO^'^K and initial value of dark energy density ^Acw = lO^GeV^. 

Moreover, it is important to mention that existence of riccions are 
possible at energy scales where higher-derivative terms in the gravita- 
tional action (7.1.1a) has significant role compared to Einstein-Hilbert 
term. At energy scales below Mew, Einstein-Hilbert term dominates 
higher-derivative terms in the action (7.1.1a). So, riccions has no di- 
rect role in the cvolutionof proposed model of the universe, but it has 
two very important contributions to the observable universe as value of 
dark energy density p^^^ — lO^GeV^ and cosmic background radiation 
temperature = 9.1 x 10^^ K at cosmic time t — 0. 
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Astronomical observations have compelling evidences that the cur- 
rent universe is dominated by dark energy (DE). The present cosmic 
dark energy density is very low, but it used to be very high in the early 
universe. The fall of dark energy density from very high to extremely 
low value can be explained if it is time-dependent. So, like other other 
cosmic dark energy models, here also, dark energy density pA is slowly 
varying function of time. 

In 1933, Bronstein proposed that, in the expanding universe, dark 
energy density decreases due its decay as a result of emission of dark 
matter or radiation [70, 71]. The radiation, so emitted, could disturb 
spectrum of 3K-microwave background radiation. So, Bronstein's orig- 
inal idea was modified and it was introduced that dark energy (DE) 
could decay to hot or cold dark matter without any harm to spectrum 
of 3K-microwave background radiation [70] . Following this idea,here it 
is demonstrated that pA{t) decays to hot dark matter(HDM) till tem- 
perature is high and cold dark matter after decoupling of matter from 
radiation. As a result, pA{t) falls from high value lO^GeV^ to currently 
low value 0.73pcr,o (where pcr,o is the critical density). This approach 
also provides a solution to cosmic coincidence problem. 

Here, it is proposed that our 4-dimensional dynamical universe grows 
below Mew and begins with topology having the distance function 

dS"^ = dt^ - (t) [dx"^ + dy^ + dz'^] (7.4. 1) 

for the spatially flat model of the universe supported by recent experi- 
ments [56, 57, 58, 81]. This space-time is a special case of hypersurface 
of the line element (7.1.1a). Here a{t) is the scale factor. 
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It is shown above that , in the beginning, the universe was very hot 
due to the background radiation (released during phase transition). 
Radiation energy density falls as 

_ Pew(r)Qew (7 4 2) 

with growing scale factor. 

Matter remains in thermal equillibrium with the background radia- 
tion for sufficiently long time. According to WMAP [81], decoupling 
of matter from background radiation takes place at td ~ 386kyr = 
1.85 X lO^^GeV"^. Equation of state for radiation is a; = 1/3. So, here, 
it is proposed that dark energy decays to hot dark matter (HDM) till 
decoupling time obeying cu^ = 1/3. 

When t > td, it decays more to cold dark matter (CDM) ( which is 
non-baryonic and pressurcless) with ujm = as well as HDM. So, ratio 
of densities of HDM and CDM is extremely small below ta- 
in what follows, development of the universe is probed using these 
ideas. 

(a) Decay of dark energy to dark matter 
The conservation equations 

TIa^j, + Ti^j, = (7.4.3) 

yield coupled equations 

PA + 3H{1 + uja)pa = -Q{t) (7.4.4a) 

and 
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Pdm + SH{1 + LUdmjpdm = Q{t) (7.4.46) 

with pAit) and Pdm(^), being dark energy density and dark matter 
density respectively at cosmic time t. Here H = a/a,ujA = Pa/pa = 
—0.4 and cudm = Pdm/pdm • Q{t) is the loss (gain) term for DE(DM) 
respectively. 

Batchelor [76] has pointed out that dissipation is natural for all ma- 
terial fluids, barring superfiuids . So, presence of a dissipative term for 
dark matter DM is reasonable. In eq. (7.4.4b), this term is obtained by 
taking 

Q{t) = SnHpdrn (7.4.5a) 

(with n being the real number) as in [77] without any harm to physics. 
With this setting for time-dependent arbitrary function Q{t), eq.(7.4.4b)is 
obtained as 

Pdm + SH{1 -n + w)pdm = 0. (7.4.56) 

Here — npdm given by Q{t) acts as dissipative pressure for DM. As Q{t) 
is proportional to pdm, this setting does not disturb the perfect fluid 
structure as required by the 'cosmological principle'. Thus eq. (7.4.5b) 
yields the effective pressure for DM as 

P{dm,eff) = {-n + W)pdm- (7.4. 5c) 

Now connecting eqs. (7.4.4b) and (7.4.5) and integrating, it is ob- 
tained that 

Pdm = Aa=^("-'^dm-i) (7.4.6a) 
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with 

A = 0.23per,oao'^""'^ (7.4.66) 

using current value of dark matter density Pdm,o = 0.23pcr,o and oq = 
a{to). 

Connecting eqs. (7.4.4a), (7.4.5) and (7.4.6a) and integrating, it is 
obtained that 

Pa = 3^+^,) , Pdm (7.4.7a) 

where 

D = 10^a34i+^A)GgV^ (7.4.76) 
using the initial condition pdm = at i = 0. 
(b) Expansion of the universe, when p\ > pdm 
Priedmann equation is given as 

/d\2 SttGn, . StcGn , . /^T . ON 

\a) ^^—(P^ + P^^ + Pr)-^—{pA + Pdra) (7.4.8) 

as (pa + Pdm) dominates over pr, which is clear from eqs. (7. 4. 2), (7. 4. 6) 
and (7.4.7). 

In this case, the equation (7.4.8) reduces to 



\a) ~ 3a3(i+-A) ^' 



yielding the solution 



a{t) = a^ 



,MGnD t 



3 9/10 



1 + 2.89 X IQ-^H 



10/9 



(7.4.10a) 
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with a;A = —0.4, Ocw = = 0) and 



SttGnD 



2.89 X lO-^^aia^GeY 



(7.4.106) 



( Gn = Mp^, Mp = lO^^GeV). 

It shows an accelerated expansion of the universe from the beginning 
itself as a/a > 0. 

Connecting eqs. (7.4.7) and (7.4.10) 



Pa 



10« 



n 



1 + 



SttGj 



n + a;A — a;. 



(7.4.11) 



But, the present universe obeys the condition 



^^A,o + f^m,o = 1, (7.4.12a) 

where Oni,o = ^^dm,o + ^^r,o- Here VLq = po/Pcr,o and pcr,o = SH^/SttGn 
{Hq being Hubble's constant) for the present universe. With these 
values the present critical density is calculated as 



p,, o = 3H^/8nGN ^ 1-2 x 10-^^GeV^ (7.4.126) 

using Hq — ho /to with ho — 0.68 and the present age of the universe 
to = 13.7Gyr = 6.6 x lO^^GeV"^ [81]. 

As, in the present universe, pm is dominated by cold dark matter 
density, cUm = is taken in eq.(6.11). Now, conecting eqs. (7.4. 11) and 
(7.4.12), it is obtained that 



n = 0.47 (7.4.13) 

Using eq.(7.4.10), the scale factor a^ at the decoupling time t^ — 
1.85 X lO^^GeV"^ and the present scale factor ao are obtained as 
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ad ~ 1.34 X lO^^Oew (7.4.14) 

and 

oo ~ 1.58 X lO^Vw (7.4.15) 
using the present age of the universe given above. 

The equation(7.4.10) shows an accelerated growth of the scale fac- 
tor when Pa > pdm exhibiting non-adiabatic expansion of the universe. 
It implies non-conservation of entropy of the universe, which is in con- 
trast to the decelerated adiabatic expansion of SMU, when cosmological 
model is radiation-dominated or matter-dominated. 

In what follows, scale factor dependence of temperature and en- 
tropy is obtained in an emperical manner, based on current temper- 
ature of the microwave background radiation Tq = 2.73K, the ini- 
tial temperature Te„ = 78.5GeV = 9.1 x lO^'^K given by eq.(5.9) and 
Qq ~ 1.58 X lO^^Cew given by eq.(7.4.15). Thus , the required emperical 
relation is obtained as 

r T 103/50 f, 

T/tJ =^ (7.4.16) 
L J a 

which yields the decoupling temperature as 

ra -1 50/103 

Td = — = 1740.42X, (7.4.17) 

L ttd J 

using eqs. (7.3.7) and (7.4.14). This value is much lower than Td = 
3000iir, obtained in the standard big - bang cosmology. These drastic 
changes are due to dominance of the dark energy in the proposed 
model. 

Using eq. (7.3.6b), entropy of the universe is calculated as 
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a^T\ (7.4.18) 

90 ^ ^ 

Current value of entropy 5*0 is supposed to be 10^^. "How could so 
high entropy of the universe be generated ?" is an old question. A so- 
lution to this problem was suggested in a seminal paper on inflationary 
model of the early universe by Guth [82] and modified version of the 
same by Linde [83] and Albrecht and Steinhardt [84]. Here, a differ- 
ent answer to this question is provided on the basis of results obtained 
above. 

Connecting eq.(6.15) and eq.(6.18) and current vlues of entropy as 
well as temperature Tq — 2.73K of the universe, it is obtained that 



aew = 0.25 (7.4.19) 

as 

77r2 

^0 = 10^7 = _(1.65 X 2.73 X lO'^aew)' 
90 ^ ^ 

Connecting eqs.(6.16) and (7.4.18) it is obtained that entropy grows 
with the scale factor a{t) as 



S = ^-^T^y^T^'a'"'/''' (7.4.20a) 

with the initial value 



,5ew = 9 X 10^^ (7.4.206) 
The equation (7.4.6) yields the rate of production of HDM 

1 59 

Phdm = -5.9 X 10-6^^ (7.4.21) 
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using ujdia =1/3 for hot dark matter. This equation shows that rate of 
production of HDM increases with growing scale factor. So, production 
of HDM is responsible for increasing entropy in the universe. In SMU, 
number of photons decide entropy of the universe. But, in PM, photons 
and HDM both are responsible for entropy. As energy of HDM increases 
due to its production, owing to decay of DE, entropy is generated in 
this model. It is unlike SMU, where entropy remains conserved. 

Connecting eqs. (7.4.6) and (7.4.16), temperature dependence of pdm 
is obtained as 

3(n-Wdm-l) 

— n no 

Putting a;dm = 1/3(0) for HDM (CDM), in eq. (7.4.22), ratio of CDM 

and HDM densities, pcdm and phdm, is obtained as 



Tew/T 



(7.4.22) 



Pcdm/ Phdm ^ew -^w/-^ 



2.06 



(7.4.23) 



using n = 0.42 from eq.(7.4.13). Eq. (7.4.23) shows that creation of 
HDM is higher at high temperature. But as temperature falls down, 
creation of CDM supersedes the production of HDM ( production of 
dark matter owing to decay of dark energy) . For the current universe, 
the ratio of eq. (7.4.23) is obtained as 



Pcdm,o/Phdm,0 = 2.1 X lO'^ (7.4.24) 

using numerical values of aewjT'ew and Tq (given above) in eq. (7.4.23). 
It shows that currently, Phdm is almost negligible compared to Pcdm- 

The scale factor Ogd , upto which dark energy vanishes due to decay, 
are obtained as 
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Osd = 6.78ao (7.4.25) 

connecting eqs.(7.4.6),(7.4.7) and (7.4.10) as well as using pA = 0. Here 

t^dm = is taken as for t > to, dark matter content is expected to be 
dominated by CDM. 

The ratio of densities of dark energy, p\, and dark matter, pdm, is 
obtained as 

47 

PA/Pdm = 10(ao/a)°-2i - —. (7.4.26) 

It is interesting to note, from eq. (7.4.26), that the gap between 
dark energy density and dark matter density decreases with the grow- 
ing scale factor. Connecting cqs. (7.4.6) and (7.4.7) and using n from 
eq.(7.4.13). This equation shows that Pdm/ Pa < 1 for 

aew < ot < 3.44ao and 0<t< SMto. (7.4.27a, b) 

The result, given by eq. (7.4.27b), solves the cosmic coincidence problem 
[24] . This approach, for solution of coincidence problem has an advan- 
tange that no scalar field is required to represent the dark energy, 
rather it is the contribution of higher-dimensional higher-derivative 
gravity to the observable universe. It is unlike the case of work in 
[75, 77], where assumed quintessence scalars are used. 

(c) The universe, in case pA < pdm 

Like eq. (7.4.27), employing the same procedure, it is also possible to 
find that Pdm/ Pa > 1 for 



a > 3.44ao and t > 3.03^0- 



(7.4.28a, b) 
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So, for t > 3.03to, the the Friedmann equation is written as 
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(-] = ^^(PA + Pdm + Pr) ^ — 3^Pdm ^ 0.27H^{ao/ af '^ 

(7.4.29) 

using eqs.(7.4.6),(7.4.8), (7.4.13) and definition of pcr,o (given above). 
Eq.(6.29) yields the solution 

r cT/inn 1 100/87 

<t) = [4<r + 0.43 X 10-^'{t - tA<dm)J , (7.4.30) 

where aA<dm — 3.44ao and iA<dm — 3.03io- Eq. (7.4.30) also shows an 
accelerated growth of the scale factor, even though p\ < pdm- It is 
interesting to see that this expansion is faster than the expansion in 
the interval < t < 3.03to- 

From cq. (7.4.25), the scale factor upto which dark energy vanishes 
is — 6.78ao. So, eq. (7.4.30) yields the corresponding time as 



4d = 3.7^0 = 50.69Gyr. (7.4.31) 

It shows that decay of dark energy will continue in the interval 
3.03^0 < t < 3.7to also. So, during this time interval, the dark matter 
will follow the rule, given by eq. (7.4.6). 

Using cqs. (7.4. 16), (7.4. 18), (7.4.25)and (7.4.31), T^d and S^a are cal- 
culated as 

/ an \ 50/103 

Tsd = To I — ) = 0.395To = 1.08K. (7.4.32a) 

Vasd/ 

Dark matter density at t = tgd is obtained as 



Pdm(sd) = 1.32 X lO-^^GeV^ 



(7.4.33) 
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using ao and asd in eq. (7.4.6). 

As dark energy vanishes when t > tsd, content of the universe will be 
dominated by CDM ,which is pressureless non-baryonic matter obeying 
the conservation equation 

Pdm + 3-ffpdm = 0. 
This equation yields scale factor dependence of pdm as 



Pdm = Pdm(sd) (asd/ a(t)j = 



3 2.54 X 10 



39 



(7.4.34) 



for t > 3.7to. 

Beyond the age of the universe 3.7to, the Priedmann equation looks 
like 



0'- 



8nG 



N 



-Pdm — 



212.8 



(7.4.35) 



using eq. (7.4.34). 

Eq. (7.4.35) yields the solution 



a{t)= a%^±U.6{t-tsd) 



2/3 



(7.4.36) 



with Osd and tgd given by eqs.(6.25) and (7.4.31) respectively. 

On taking (+) sign in eq. (7.4.36), decelerated expansion of the uni- 
verse is obtained beyond t > 3.7to continuing for ever. But the (-) 
sign, in eq.(6.36), exhibits a contracting universe beyond t > 3.7to. 
So, ultimately the contracting universe is expected to collapse to a 
very small size with scale factor, possibly equal to = 0.25. Prom 
eq. (7.4.36), with (-) sign, the collapse time is calculated as 
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a' 



3/2 



^col — ^sd ~l~ 



14.6 



- 18.08^0 = 247.73Gyr 



(7.4.37) 



using Osd = 6.78ao and tsd — 3.7io- 

7.5. Elementeiry peirticles, Primordial nucleosynthesis and 

Structure formation 

(i) Creation of particles 

In SMU, it is assumed that elementary particles such as leptons, 
mesons, nucleons and their anti-particles were produced at the epoch 

of big-bang. Here, production of scalar and spin-1/2 particles is pro- 
posed in the very early universe due to topological changes caused by 
expansion of the universe. 

(a) Creation of spinless particles 

The scalar field 4> obeys the Klein-Gordon equation 



where is mass of 0. Expanding in terms of mode k, eq. (7.5.1) is 
written as 



(□ + m!)0 = O, 



(7.5.1) 



(1 + 2.89 X 10'^H)(j)k + 



28.9 



X 10-^^0+ mi(l + 2.89 X lO^^^t) 



3 



_^(l + 2.89xl0-i6t)-"/'^ 

using the scale factor a{t), given by eq.(6.10a). 
For small t, eq. (7.5.2) is obtained as 




(7.5.2) 




(7.5.3) 
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r = 1 + 2.89 X IQ-^H. 
Eq.(7.5.3) yields the normahzed solution 



(7.5.4) 



[2(2.)3V6^]-V^e-^[l + i::^^^^], (7.5.5a) 



where 



and 



20(10^)^ 



(7.5.56) 



62 = 10=^'(m| + 



9^ 



Here iFi(a,6, c) is the confluent hypergeometric function. 
For large t, eq. (7.5.2) reduces to 



dr"^ 3 dr 
which integrates to 



(7.5.5c) 



(7.5.6) 



out ^ ^-7/6 



CiJ_7/6(^t) +02^-7/6 (^T 



where 



(7.5.7a) 



A = lO^^-^m^. (7.5.76) 

and Jp{x) and l^(x) are Bessel's function of first and second kind re- 
spectively. For large x, Bessel's functions can be approximated as 

cos{x — 7r/4 — p7r/2) 



Jp{x) ~ 



y/Trx/2 



(7.5.7c) 
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and 

Using these approximations, when r is large, eq. (7.5.7a) looks like 

.-5/3 



iOUt 



-fe - -j===[cicos{Ar + -) + C2sin{Ar + -)] (7.5.8) 

Solutions (7.5.5a) and (7.5.8) yield number of produced spinless par- 
ticles (for mode k) per unit volume as 

-10/3 

W = ^j^\X\'^0, (7.5.9a) 



where 
X = 



^ _ 3&i\^r _ 2 cos(At) - A-y^rsin(AT) + i | y^r 

-^n('-^)^-H^^)-(^-w)-^-"H- 

(7.5.96) 

is defined in Appendix B. Non-zero shows creation of spinless 
particles. 

(b) Creation of spin- 1/2 particles 
The spin- 1/2 field ip satisfies the Dirac equation 



(lYD^-mfyiP = 0, (7.5.10) 

where m/ is mass oiip, are Dirac matrices in curved space-time and 
are covariant derivatives defind above, ip can be written as 

^ = XI ^iKs'^ik,s + dl^^ipnk,s) (7.5.11) 

s=±l k 

with hk,s and dk,s-i given in Appendix B. Eqs.(7.5.10) and (7.5.11) yield 
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(7.5.12a) 



iiYD^-mf)ijjrk,s^O. (7.5.126) 

Now using the operator {— i'-y^ 0^ — 171 f) from left of eqs.(7.5.12a) and 
(7.5.12b), it is obtained that 



;n + -i? + mj)vi = o, 



where i) = ipik,s{i^iik,s)- 
Writing 

and 



■ipiik,s = fiik,s{t)e'''-^Us, 
eq.(7.5.13) is obtained as 



(7.5.13) 



S~f , 10 , 



+ 



+ lO'^ 



k'^ 4.93 X 10 



-15 



11/9 



dr^ 3 dT 
where f ^ fik,s{fiik,s)- 

For smaU t , eq.(7.5.14) looks like 



3t 



/=0, (7.5.14) 



d^f 10 df 



^ 2 _|_ ■ 



r — 



20^ 

Qfl2 



/ = 0. (7.5.15) 



This equation is like the equation (7.5.4a). So, its solution has the 
same form. Using this solution in ipik,s defined above , it is obtained 
that 



!T\,'h2 



a/2(2^0V^ 
For large t , eq.(7.5.14) reduces to 



^ , (-?A + 10/v^) 

i H ' 

10 



(7.5.16) 
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d^f 10 df 



+ lO^^mJr/ = 



(7.5.17) 



yielding 



^out 



-5/3 



IIk,s 



■cos(/lr)e*''-^Ms. 



(7.5.18) 



Using ipYk,s ill the definition of Pk,s, given in the Appendix 



B, 



= (27r)V-V3 



cos{At) e 



■ , (-361 + 10^75^) 

^ + 10 \ 



(7.5.19) 

So, the number of created spin- 1/2 particles per unit volume is ob- 
tained from eq.(7.5.19)as 



\Pk,s\' - 



(2vr) 



zT 



-10/3 



V , h 2 



(7.5.20) 



2^^^A^/b2 

Eqs.(7.5.9a,b) and (7.5.20) show creation of spinless and spin-1/2 
particles due tochanging gravitational field in the proposed speeded-up 
model. The scale factor, given by eq.(7.4.10a) show that the significant 
change in the gravitational field is possible , in this model, when 



t > 3.46 X lO^^GeV"^ - 7.2 x lO'^sec. (7.5.21) 

as a{t) remains almost constant upto this epoch. So, here, particle 
production is expected, when the universe is around 2.3 x 10~^sec. old. 
Also, these results show that particle production falls down rapidly as 
time increases. Thus, like other models, here also particle production 
is expected in the early stages of the universe. 

7.2 Primordial Nucleosynthesis 
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Hydrogen (H) is the major component of baryonic matter in the 
universe. The next main component is Hehum-4 {^He) . Occurrence 
of other hght elements and metals is very small. It is found unlikely 
that abundance of ^ife, Deuterium {D) and other light elements, in 
the universe, being caused by burn out of H in stars [85, 86] So, it is 
argued that the required amount of these might have produced in the 
early universe. 

At the epoch of its formation, helium production depends upon neu- 
tron {n) concentration, which is determined by weak interaction reac- 
tions given as 

n + u^p-\-e~,n-\-e^^p + i' (7.5.22) 

where p stands for proton and v for neutrino. This chemical equillib- 
rium is maintained till weak reaction rate >> H, where H is the 
expansion rate of the universe given by cq. (7.4.9) and r„ ~ 1.3G|nT^ 
with Fermi constant Gp = ira^/V^M^ = 1.17 x IQ-^GeV^ With the 
expansion of the universe, temperature decreases so weak interaction 
rate slows down . As a result, at the freeze-out temperature T*, 

Fw ~ H. (7.5.23a) 

Connecting eqs.(7.4.7b), (7.4.9), (7.4.16) and (7.5.23a), is ob- 
tained that 

1.3 X (1.17)2 X IQ-^^T^ ~ 2.89 x 10"^^ (^^^'^ 
which yields 

= 0.9MeV. (7.5.236) 

Eqs.(7.4.10a,b) and (7.4.16) yield time dependence of temperature , 
for this model, as 
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-'MeV 

Using this result, the freeze-out time t^, corresponding to T*, is ob- 
tained as 

U ^ 5.55sec. (7.5.25) 

In SMU,r* ~ 0.86MeV and U ^ Isec [87]. It is obtained that, in 
PM, is a bit higher and freeze-out takes place later. The reason 
for these differences is the basic difference in developement of these 
models. SMU is driven by elementary particles present in the early 
universe and it expands adiabatically. The proposed model is driven 
by dark energy, with very high density at the beginning and it expands 
with acceleration. 

Like [88], in this subsection, temperature Tg is measured in units 
lO^K, so IMeV becomes 11.6 in Tg temperature. Now, eq.(7.5.24) 
looks as 



474 

(7.5.26) 



For large temperature Tg > 10, neutron abundance is given by 



1 + e 



-1 



(7.5.27a) 



where Qg = m„ — rUp — 15. At freeze-out temperature ^ 0.9MeV, 
which is equivalent to Tg = 10.44, 



X* = 0.19. (7.5.276) 

Due to slightly higher freeze-out temperature, in the present model, 
neutron concentration is obtained higher than X* = 0.16 in SMU. 
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When Fw << H i.e. when chemical equiUibrium between n and p 
freez-out, neutron concentration is determined through neutron-decay 

n^p + e~ + P. (7.5.28) 

Thus, for t > t^, neutron abundance is given as 

Xn(t) = 0.19e-*/^, (7.5.29) 

where r = (885.7 ± 0.8)sec. is the neutron hfe-time. 
First step, in the formation of complex nuclei, is 

p + n D + J 

showing chemical equiUibrium between deutron (D) , nucleons and 
photon (7). This equiUibrium gives deutron abundance by weight [86] 
as 

Xd = 0.2 X W-'\nBh^)Xr,XpT^/^exp{Bn/n), (7.5.30) 
where Xd = 2710/^5, X„ = nn/nB,Xp = np/nB,Xn + Xp ^ 1, deu- 
terium binding energy Be> = mp + m„ — ttid — 2.23MeV which yields 
Bog = 25.82 (in units of Tg) and [VlBh^) is the baryon number density. 
Light heavy elements '^He is formed through reactions 

n + p^ D + -f, (7.5.31a) 

{i)D + D He + n, {ii)D + D H + p (7.5.316, c) 

and 

{ii)D +3 H He + n, {lifHe + D He + p, (7.5.31d, e) 
where stands for tritium. 
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These reactions show that sufficient deuterium abundance is required 
for nucleosynthesis of '^He, which is supposed to be around 25% of the 
baryonic content of the universe. But after weak-interaction freeze out, 
deuterium abundance is not sufficient unless temperature is very low. 
For example, at Tg = 5.8(r = 0.5MeV),XD = 2 x lO^^^ is obtained 
using eq.(7.5.30). 

Rates for reactions (7.5. 31b, c) are given as [88, Appendix] 

ADD(i) = 3.97 X lO^Tg^^/^e-^-^^^/^s^' [1 + 0.098r9 + 0.876^9 
+O.6T9 - 0.041T9/^ - 0.071X9/^] (f^B/i=^)sec-^ 

ADD(n) = 4.17 X 10^r9"^/V^-^^^/^9^'[l + 0.09879/^ + 0.518^9/^ 

-0.355^9 - 0.01079"^/^ - 0.018X9/^] (i7B/i2)sec-^ 

(7.5.32a, 6) 

Making approximations in eqs.(7.5.31a,b), it is obtained that 

^DD — ^DD{i) + ^DD{ii) 

~ 8.14 X 10'T;'/'e-'-'''/''^'\nBh')sec-' 

(7.5.33) 

Mukhanov [87] has obtained the condition for conversion of sufficient 
amount of D to ^He and as 

^XDoXot ~ 1 (7.5.34) 
From eqs.(7.26), (7.29)-(7.32), it is obtained that 

I 25.82 4.258 | 

6.43 X IQ-^^-^eV ^« ^^^7 (l^s/i^)^ ~ 1 (7.5.35) 

showing a relation between Qsh^ and Tgj. Tgi is the temperature at 
which nucleosynthesis takes place. For Q^/i^ = 0.05, this equation 
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yields 

T^i ~ 1.16, (7.5.36) 

which corresponds to Tj — O.lMeV. It shows that if baryon number 

density fi^/i^ is ~ 0.05, nucleosynthesis oi'^He may take place at tem- 
perature Tj = O.lMeV. In SMU, this temperature is ~ 0.086MeV. 
Eqs.(7.29) and (7.35) yield deuterium concencentration at Tj as 

Xoi ^ 8.93 X 10"^ (7.5.37) 

It happens at the cosmic time 

U = 357.5sec. (7.5.38) 

Thus, it is obtained that, in the present model, nucleosynthesis of ^He 
begins much later compared to SMU, where it happens at time ~ 
lOOsec [87]. But, it is still less than life-time of a neutron, which is 
required for the necleosynthesis. 

The final ^He abundance is determined by available free neutrons at 
cosmic time ti. As total weight of ^iJe is due to neutron and proton, 
so its final abundance by weight is given by [88] as 

Xi = 2Xlexp{-tf It) = 2Xlexp{-Q.b?,b/ {Tif-') (7.5.39) 

using eqs. (7.5.26) and r — 886.5sec. Here a stands for ^He. 
The rate of production of ^ife is given by the equation 

= 2\DDin)Xl. (7.5.40) 

As free neutrons are captured in ^iJe, helium production dominates 
neutron decay. It happens when [88] 



T 



(7.5.41) 
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Approximating \DD{n), given by eq.(7.5.32b), around Tg = 
obtained that 
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1.16, it is 



ADD(n) ^ 6.98 X lO^Tg ^^^exp{-4:.258/Tg^^){nBh'^). (7.5.42) 
Eqs. (7.5.30), (7.5.39)-(7.5.41) yield 



— ~ 6.98 X 10 To' e ^9 x 0.19e ^9 
x(0.2 x 0.81 X lO'^YinBh^f, 

which implies that 



Tf = ^ll^i (7 5 43) 

' 37.96 - l/nTg - 3ln{nBh'') + iff + ^'■^■^"'^ 

At T = Tgi = 1.16, 



J ^ 51.64 

' 42.08 - 3/n(l]B/i2)" I ■ ■ J 

Using eq. (7.5.44) in eq. (7.5.39), helium abundance by weight, as 
function of baryon number density {flBh^)M obtained as 

f r r 42.08 - 3Zn((]B/i2)-ii.9i 
~ 0.38exp[- 0.535| 5164"^ | J" (7.5.45) 

For Q^/i^ = 0.05, it is obtained that 

Xl ~ 0.23. (7.5.46) 

For T < T/, neutron abundance is governed by X„ = -2X« with 
Xa given by eq. (7.5.39). So, 

— —'^^DD{n)Xl), 
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where Xd and \DD(ii) ai'e given by eqs. (7.5.30) and (7.5.41) respec- 
tively. Using eq. (7.5.26) and integrating, this equation yields neutron 
abundance below T/ as 



(7.5.47a) 

which shows that concentration of free neutrons drops very fast. It is 
due to capture of free neutrons in ^ife. Prom eq.(7.5.30), it is obtained 
that neutron and deuterium concentrations become equal at tempera- 
ture Tg given as 



' 29.45 - Ininsh'^) - UnTf {'-^-^'O) 



Deutron concentration is given by [88] 



Xn ^ -2XDDi^^)Xl = -2 X 8.14 X 10'T-'/'{QBh')e''-'''/<'xl, 
which integrates to 



X 



D 



'(l/Xn{T,^)) + (f)B/i')e'°-^^{(T9'^)-'-'e-^-2^«/(^9)''' 

-(r9)-2-57e-4-258/(T,)V3^1-l^ 

(7.5.47c) 

where XD{Ti)=X^{Ti). 

Eqs. (7.5.47a) and (7.5.47c) show that, below T/, abundance of free 
neutron and deuterium drop rapidly in PM also like SMU. 
(iii) Growth of inhomogeneities in the proposed model 

Por small inhomogeneities, linear perturbation of Einstein equations 
is enough. In this case, contrast density 5 — 5pdm/pdm (with pdm, being 
the cold dark matter energy density and 5pdm is small fluctuation in 
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Pdm presenting inhomogeneity) , obeys a linear differential equation in 
the homogeneous model of the universe [89]. Due to linearity, 6 can be 
expanded in modes k as 



(7.5.48) 



For modes k with proper wavelength F < H ^(t){H ^(t) being the 
Bubble's radius), the perturbation equation looks like 



Ok + 2-dk + 



- A7rGpam)Sk = 0, (7.5.49) 



a \ 

which is free from gauge ambiguities [89, eq.(4.158)]. So, Sk depends 
on t only. 

Eqs. (7.4.5c) and (7.4.13) give effective pressure for DM as 



It yields 



P(dm,efr) = (-0.47 + Wdm)Pdm- 



2 C^P(dm,eff) n 4^ , 

V, — — ^ — —0.47 + w. 



s 



dp, 



'dm 



(7.5.50) 



(7.5.51) 



Connecting eqs.(7.4.6), (7.4.13), (7.5.49) and (7.5.51), the differential 
equation for 6k is obtained as 



(-0.47 + Wdm)^ - 0.345i/o'aJ-'7o'^°-^^+^'^-^ 



5k + 2-5k + 



with pcr,o = 3Hq/87iG. 

Using eq.(7.4.10), eq.(7.5.52) looks like 



4 = 

(7.5.52) 



(PSk , 20 dSk , ^„3i 



0.345//2aJ-59 



-0.47 + Wdm)- 



.11/9 



^10(0.23+Wdni)/3_ 

where r is defined in eq.(7.5.4). 



0, 



(7.5.53) 
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As mentioned above, for CDM Wdm = 0, so eq. (7.5.51) is approxi- 
mated as 

<P5, 2Qd5k 3.45 X W^Hlal^^ ^ 

+ rm^^ — = 0' (^■^■^4) 

which is integrated to 



^-11/18 



(Cl + C2) Jll0/2l(2i6T^/^°) + i(ci - C2)yil0/2l(2i6T^/^°) 

(7.5.55) 

where i = and 6^ = 3.45 x IQ^^Hlal^^. Here //^^^(a;) and h!^^\x) 
are Hankel's functions as well as Jp{x) and l^(x) are Bessel's functions. 
Moreover c\ and C2 are integration constants. 

As discussed in section 7.4, CDM dominates HDM in the late uni- 
verse, so structure formation is expected for large r. Using approxima- 
tions (7.7c,d)for Bessel functions, for large r in eq.(7.5.55)and doing 
some manipulations, it is obtained that 

5,f«r-°-6V^''", (7.5.56) 
which shows growth of structure formation as cosmic time t increases. 

8 Summciry of Results 

Thus it is found above that dual role of the Ricci scalar as a geometri- 
cal field as well as a scalar physical field with its particle concept termed 
as riccion (which is different from 'scalaron' or 'graviton') reflects many 
important features of gravity and cosmology. Interestingly, these fea- 
tures are inspired by geometry of the space-time. In section 3, it is ob- 
tained that riccion behaves like an instanton. The instanton solution. 



DUAL NATURE OF THE RICCI SCALAR 107 

SO obtained (using dual role of the Ricci scalar), causes inflation prior to 
GUT phase transition. This inflationary scenario is called 'primordial 
inflation'. Riccions, emerging from higher dimensional space-time, de- 
compose into fermionic and anti-fermionic counterparts, called riccions 
and anti — riccions when parity is violated. As a riccion and a riccino 
have the same mass, these arc supers- symmetric partners too. Mass of 
riccions and riccinos, being inversely proportional to gravitational con- 
stant, is very high. So, during the course of expansion of the universe, 
these might have decayed due to emission of observed elementary par- 
ticles. If extra-dimension D of the (4 -|- £)) —dimensional space-time is 
energy mass scale dependent, results of one-loop rcnormalization yield 
D as positive real number above the energy scale 1.76 x lO^^GeV, not as 
positve integer. It supports fractal dimension, not topological dimen- 
sion above this scale. But when D is taken independent of energy mass 
scale, renormalization group equations for riccion yield D — 6, making 
dimension of the space-time -|- 4 = 10 as required in super-string 
theory. 

Riccion yields many interseting cosmological consequences too. Ho- 
mogeneous and inhomogeneous models are dervied here using dual role 
of the Ricci scalar. It is shown that the universe bounces at the critical 
temperature ~ 10^'^GcV supporting a non-singular biginning of cos- 
mology. Creation of spinless and spin-1/2 particles in these models are 
discussed here. Lastly, a new cosmological scenario, inspired by dual 
nature of the Ricci scalar, is proposed here. It is discussed at length in 
the preceding section and briefly it is as follows. 

In contrast to the original big — bang theory, the proposed cosmol- 
ogy answers many basic questions as (i) "What is the fireballT , (ii) 
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"How docs it burst out?'', (iii) "What is the background radiation?" 
and "What are initial values of temperature and energy density?" 
Moreover, the proposed cosmological model is free from initial sin- 
gularity. It is found that the dark energy violates the strong energy 
condition showing "bounce" of the universe, which is consistent with 
singularity — free model of cosmology [16, 21, 22]. The initial scale 
factor is computed to be a^^ = 0.25. It is in contrast to the standard 
model of cosmology SMU, which encounters with singularity having 
zero scale factor, infinite energy density and infinite temperature. 

The present value of dark energy density is supposed to be ~ 7.3 x 
10^'^^Gev*^, which is 53 orders below its initial values lO^Gev"^ . The 
question "How docs dark energy falls by 53 orders in the current uni- 
verse?" is answered by the result, in section 6(a), showing that dark 
energy decays to dark matter, obeying the rule, given by eq.(6.7). 
Upto the decoupling time td — 386kyr, matter remains in thermal 
equilibrium with the background radiation, so produced dark matter 
upto td is supposed to be HDM. But when t > td, production of CDM 
is more than HDM. The current ratio of HDM density and CDM den- 
sity is found to be 5 X 10~^°. In section 6, it is found that creation of 
HDM raises entropy of the universe upto ~ 10^^. During the phase of 
accelerated expansion, temperature falls as a^t)^^'^^^*^^ . 

It is found that one of the two types of dynamical changes of the 
universe are possible, beyond 3.7to — 50.69Gyrs, (i) decelerated expan- 
sion and (ii) contraction. If future universe expands with deceleration 
, it will expand for ever. But in the case of contraction, it will collapse 
by 247.73Gyrs. 
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As dark energy dominates over dark matter^ from the beginning 
upto the time 3.03to, no cosmic coincidence problem arises in the pro- 
posed scenario. Moreover, in the preceding section, some other basic 
problems such as creation of particles in the early universe, primordial 
nucleosynthesis and structure formation in the late universe. In the 
proposed model, it is shown that particles are created due to topologi- 
cal changes which is unlike the SMU, where it is assumed that elemen- 
tary particles are created at the time of big-bang. It is also shown, in 
section 7.2, that process of primordial nucleosynthesis goes very well in 
the proposed model predicting ~ 23% Helium-4 abundance by weight. 
In section 7.3, it is shown that ,in the late universe, inhomogeneities 
in CDM grow exponentially with time. Thus the proposed model is 
able to provide possible solutions to many cosmological problems with 
prediction for the future universe. 

Thus, it is found that dual nature of the Ricci scalar has enough 
potential to reveal many interesting aspects of gravitation and cosmol- 
ogy. It many yield some interesting results, if riccion is obtained from 
brane-gravity and Einstein - Gauss- Bonnet theory. 

Appendix A 
Riccion and Graviton 

From the action 




d'^xd^y^-gi^i+D) 



^(2+D)D 



IGtt 



+ «(4+£»)-R(4+D) + 



7(4+D)(-R(4+D) - 




(D+4)-R(D+4)) ) 



(Al) 



the gravitational equations are obtained as 
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M(2+-D) 1 (2) 

{RmN --^gMNR{4+D))+Ol(4+D) Hj^ jy + 7(4+D) -^M Af = 0' (A.2a) 

where 



{A.2b) 

and 

(2) 2 2 6(i^ + 3) , 1 2 

-f^MAT = 3it!.jv^^-3cj'MArn(4+D)-R(4+D) " (D - 2) ^~2^^'^^^'^+^^ i^+D) 

+2n(4+D)-R(D+4)-RMA?" + -RfMAr} ~ Q^dMNRf^+D) + '^R^{i-\-D)RMN ■ {A.2c) 

Taking qmn — Vmn + hMN with rjMN being (4 + D) -dimensional 
Minkowskian metric tensor components and Hmn as small fluctuations, 
the equation for graviton are obtained as 



□(4+D)^MiV = {A.3) 

neglecting higher-orders of h. 

On compactification of M'^ (g) to M^, eq.(A.3) reduces to the 
equation for 4-dimensional graviton as 

for the space time 



^5-2 = g^^dz'^dx" - p^[del + sin^OidOl H h sm^6'i ■ ■ ■ sin'^9(D-i)d9l). 

{A.5) 

The 4-dimensional graviton equation (A. 4) is like usual 4-dimensional 
graviton equation ( the equation derived from 4-dimensional action) 
only for Z = 0. Thus the massless graviton is obtained for Z = only. 
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As explained, in section 2, the trace of equations (A. 2) leads to the 
riccion equation 



[n + ^^R + ml + ^R^]R + ^^0, (A6a) 

where 



A 

1? 



D 



2(D + 3) 

_{D+2)XVd , DRp , 1 „2 
167rG(4+D) 2(D + 3) 2^^ 
1 



4(D + 3)a' 

+ 2)AM(2+^)V^s 



-?7 



+ ^^^^ + Ir^^XRl 



167r ■ 4(D + 3) 6 

(yl.66, c, (i, e) 

The graviton /i^,^ has 5 degrees of freedom ( 2 spin- 2 graviton, 2 
spin-1 gravi- vector (gravi-photon) and 1 scalar). The scalar mode / 
satisfies the equation 



(A7) 



Comparison of eqs.(A.6) and (A. 7) show many differences between the 
scalar mode f of graviton and the riccion (R) e.g. ^, A and 'd are 
vanishing for /, but these are non- vanishing for R. Eq.(A.7) shows 
{massY for / as 



l{l + D 



(A 



whereas {massY for R (?77.^), given by eq.(A.6c), depends on G'(4+£i), Vd 
and Rd (given in section 2). = for Z = 0, but can vanish only 
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when gravity is probed upto ~ lO'^^cm. As mentioned above, gravity 
could be probed upto minimum distance 1cm only. Thus, {massy of 
riccion can not vanish. 

Thus, even though / and R are scalars arising from gravity, both are 
different. Riccion can not arise without higher-derivative terms in the 
action, but graviton is also obtained from Einstein- Hilbert action with 
lagrangian density R/SnG. 

In a locally inertial co-ordinate system 



1 d r , d \ 



showing that in the UR^ R and □ are like other scalar fields and □ 
in Klein -Gordon equation for 0, according to the principle of covari- 
ance. So, there is no harm in making loop quantum corrections and 
renormalization for R like 0. 

Appendix B 



In eq.(6.4.6), using the convolution theorem the function e**^ '' 

4/7 

can be written as 



Ak.r 



-4/7 



/•oo 
/ ' 



dye 



-ik.y 



1-^ 



-4/7 



{B.l) 



where r] has dimension of length, as used in section 2, for dimensional 
correction. Here in the integral dy introduces a quantity of length 
dimension which is compensated by 77. 
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Evaluation of the integral in cq.(B.l), is done as follows. 
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oo Q-ik.y 



poo 

/ dy 

J ro+e 



y 



4/7 



dy-, 



-ikY 



n\ r(a; + 4/7) 

^oj-3/7 ^-x{y-ro) 



[-iky 



[-iky 



2 n\ r(a; + 4/7)yo 



(ro + e)("+^/^)e'^ 



1 ui uiiui — 1) 

- H \ ^ - 

X ex^ e^x^ 



+ 



{B.2) 



neglecting terms containing (ro + e) ^ compared to terms containing 



Performing further integration in eq.(B.2), one obtains that 



dy- 



-ik.y 



ro+e 



y 



4/7 



n=0 
^3/7 



-ikY 1 



n\ r{uj + 4/7) 



(ro + e)("+^/^) X 



r(a; - 3/7) + u;r(u; - 10/7) + 
cu{uj - l)T{(jj - 17/7) + 
ikY 1 



n=0 



n\ r{u} + 4/7) 



{ro + e) 



(n+4/7)^3/7 



= (ro + e)^/V/^cosA;(ro + e), 



(S.3) 
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when uj oo. Here, symmetry under k —k caused by symmetry of 
the model under r — > — r is used. 
The same procedure yields 

r r 1 -4/7 

/ dye-'^-yy-'' 1 - ^ = (ro + eY^'^'h^/' cosk{ro + e), {BA) 

POO 

/ dye-''-yy-'/'[y - vo]-''/' = (ro + e)-'/'e-''/'cosk{ro + e), {B.5) 

J ro+e 



f oo 



PC 

/ dye'''-yy-^/''[y - rof = (ro + e)-^^h''/^cosk{ro + e), (5.6) 

J ro+e 

POO 

/ dye-"'-yy-"'/''[y-ro]-'''/' = {ro + e)-^''/'e-'^/'cosk{ro + e). {B.7) 

J ro+e 
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